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Nomenclature 


L 

®12 

“p2 


vector of pitch line 

major and minor axes of gear cradle ellipse 

major and minor axes of pinion cradle ellipse 

gear head cutter diameter 
Mean pitch cone distance 
gear ratio; mj ^2 “ 

pinion cutting ratio 

gear cutting ratio 


n 


(F) 

9 


n 


(p) 


q 

qF* qp 


unit normal to pinion and gear tooth surface 
direction angle of normal in coordinate system 
parameters of pinion and gear machine-tool setting 


r 


(P) 
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r 


(F) 


Up. Up 


(F) (P) 

V 9 V 


(1) (2) 

V > V 




radius of generating surface measured in plane X =0 
(i = 1, 2) “ 

generating cone surface coordinate 

surface E , Z contact point velocity 
r r 

surface Z 2 contact point velocity 

velocity of intersection point at gear cutter axis in plane II 
velocity of intersection point at pinion cutter axis in plane II 


W 


Yi. Y2 


'f- 


gear cutter width 
mean spiral angle 

pinion and gear pitch angle 

pinion root angle 

orientation angle of pinion 


and gear cradle ellipses 


V 
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e 

n 




aej, aLj^ 






n 


^ 1 * ^2 


S’ s 


S’ s 

4;<P) 

^C c 


( 1 ) ( 2 ) 
0) , 0) ^ 


(F) (p) 

0) 


rotation angle of frame S, relative to frame about axis Z, 

n f 1 

rotation angle of frame relative to frame S 

f n 

pinion and gear dedendum angle 
pinion machine-tool settings 

generating pinion and gear cone surface coordinates 

motion parameter of cradle ellipse 

motion parameter of pinion and gear ellipses 

plane of normals 

pinion and gear tooth surfaces 

generating surfaces of pinion and gear 
generating surfaces rotation angle 

pinion and gear blade angle 
pinion and gear angular velocities 

cradle angular velocities for cutting the pinion and gear 


Cartesian Coordinate Frame 


.(J) 

s 

,(j) 


connected to tool cone, j = F, p 


conntected to cradle 


c 

m 



connected to machine frame, i = 1, 

connected to pinion, gear 

fixed to machine, used for mesh of 
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S_ and Z- 

r i 




SUMMARY 


Spiral bevel gears are used in many applications where mechanical 
power must be transmitted between intersecting axes of drive shafts. 

Namely, two such applications are the rear axle differential gearbox for 
land vehicles and the transmissions used in helicopters. For spiral 
bevel gears, there is a continuing need for ever stronger, lighter weight, 
longer-lived and quieter running gears. Above all, a rapid and economical 
manufacturing method is essential to the Industries that use bevel gearing 
in their products. 

For many years, the Gleason Works (Dudley, 1962; Anon. Gleason Works, 
1964 and 1980) has provided the machinery for manufacture of spiral bevel 
gears. There are several Important advantages to the Gleason methods of 
manufacture over bobbing methods. The machines are rigid and produce gears 
of high quality and consistency. The cutting methods may be used for both 
milling and grinding. Grinding is especially important for producing 
hardened high quality aircraft gears. Both milling and grinding are 
possible with Gleason’s method. The velocity of the cutting wheel does not 
have to be related in any way with the machine’s generating motion. 

Generally speaking Gleason’s method for generation of spiral bevel 
gears does not provide conjugate gear tooth surfaces. This means that the 
gear ratio is not constant during the tooth engagement cycle, and, 
therefore, there are kinematical errors in the transformation of rotation 
from the driving gear to the driven gear. The research that had been 
performed by the Gleason Works was directed at the minimization of gear 
kinematical errors and the improvement of gear bearing contact by using 
special machine-tool settings. The determination of such machine-tool 
setting is accomplished by a computer program. It is known that Oerllcon 
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(Switzerland) and Klingelnberg (West Germany) have developed methods for 
generation of spiral bevel gears that can provide conjugate gear tooth 
surfaces. The disadvantage of these methods is that the gear tooth 
surfaces cannot be ground and the tooth element proportions are unfavorable 
due to the constant height of the teeth. 

The objective of the new method for generation of spiral bevel gears 
presented herein was to find a way to eliminate the kinematical errors for 
spiral bevel gears, obtain gears with higher contact ratio, and improve 
bearing contact and conditions of lubrication, while using the existing 
Gleason* s equipment and retaining all the advantages of the Gleason system 
of manufacturing such gears. The proposed method for generation is based 
on the following: 

(i) Four surfaces - two generating cone surfaces (Z and Z ) and gear 

P ^ 

and pinion tooth surfaces (Z2 Z^^) are in continuous tangency at 

every instant. The ratio of angular velocities in motion of the 
above mentioned surfaces satisfies the requirement that the 
generated pinion and gear transform rotation with zero kinematical 
errors . 

(ii) The cones have a common normal at the instantaneous point of contact 
but their surfaces interfere with each other in the neighborhood of 
contact point. 

(iii) The point of contact of the above mentioned surfaces moves in a 

plane (II) that is rigidly connected to the gear housing. The normal 
to the contacting surfaces lies in plane TI and performs a parallel 
motion in the process of meshing. 

(iv) Due to the elasticity of gear tooth surfaces their contact is 

spread over an elliptical area. The proposed method for generation 
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provides that the instantaneous contact ellipse moves along but not 
across the gear tooth surface. This provides improved conditions 
for lubrication and a higher contact ratio will result. 

(v) Until now, the reduction of kinematical errors of Gleason spiral 
bevel gears was a subject of the computer search for the optimal 
machine-tool settings. The computer program developed for this 
purpose is called the TCA (Tooth Contact Analysis) program. The 
proposed method is for direct determination of machine-tool settings 
that result in zero kinematical errors because the gear tooth 
surfaces are generated as conjugate gear tooth surfaces. 

(vi) A new TCA program directed at the simulation of bearing contact and 
the influence of errors of assembly and manufacturing has been 
developed. 

The contents of this report covers the new method for generation of 
spiral bevel gears, their geometry, the bearing contact and simulation of 
meshing. Special attention has been paid to the proposed principle of 
performance of parallel motion by two related ellipses that results in the 
desired parallel motion for the contact normal. 

1. THE GLEASON MANUFACTURING METHOD 

The gear cutter cuts a single space during a single index cycle. The 
gear cutter is mounted to the cradle of the cutting machine. The machine 
cradle with the cutter may be imagined as a crown gear that meshes with the 
gear being cut. The cradle with the mounted head cutter rotates slowly 


3 


about its axis, as does the gear which is being cut. The combined process 
generates the gear tooth surface. The cradle only rotates far enough so 
that one space is cut out and then it rapidly reverses while the workpiece 
is withdrawn from the cutter and indexed ahead in preparation for the 
cutting the next tooth. The desired cutter velocity is provided while the 
cutter spins about its axis which itself moves in a circular path. 

We consider that two generating surfaces, and are used for the 
generation of the pinion tooth surface, , and the gear tooth surface, Z 2 > 
respectively. Both sides of the gear tooth are cut simultaneously (duplex 
method) but both sides of the pinion tooth are cut separately (single 
method). The basic machine-tool settings provide that four surfaces, 

Z^, Z^ and Z^, are in contact at the main contact point. In the process of 
meshing, surfaces and Z^ , and respectively surfaces Zp and Z 2 , contact 
each other at every instant at a line (contact line) which is a spatial 
curve. The shape of the contact line and its location on the contacting 
surfaces is changed in the process of meshing. The generated pinion and 
gear tooth surfaces are in contact at a point (contact point) at every 


instants 


A head-cutter used for the gear generation is shown in Fig, 1,1. The 
shapes of. the blades of the head-cutter are straight lines which generate a 
cone while the head-cutter rotates about axis C-C. The angular velocity 
about axis C-C does not depend on the generation motion but only on the 
desired cutting velocity. Two head-cutters are used for the pinion 
generation; they are provided with one-sided blades and cut the respective 
tooth sides separately. The head cutter is mounted to the cradle of the 
machine. Fig. 1.2 shows schematically the positioning of the cradle of the 
gear machine, the head cutter and the gear to-be generated. 


i 
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2. GENERATING SURFACES AND COORDINATE STSTQIS 


The generating surface is a cone surface (Fig. 2.1). This surface is 

generated in coordinate system while the blades of the head-cutter 

rotate about axis C-C (Fig. 1.1). The generation of gear tooth surfaces is 

based on application of two tool surfaces, and Ep, which generate gears 

1 and 2, respectively. The generating surfaces (generating cones) do not 

(F) (P) 

coincide: they have different cone angles and , and different 

fF) (P) 

mean radii r^ ' and r^ ' (Fig* 2.1, a). Special machine-tool settings, AEj^ 
and ALj (Fig. 2.3,b), must be used for the generation of the pinion. 
Considering the generation of gear 2 tooth surface we use the 

(P) 

following coordinate systems: (i) which is rigidly connected to the 

( 2 ) 

generating surface (Fig. 2.1,b); (il) the fixed coordinate system 

which is rigidly connected to the frame of the cutting machine, and (lii) 

the coordinate system which is rigidly connected to gear 2 (Fig. 2.2). 

( 2 ) 

In the process of generation the generating surface rotates about the 

(P) 

- axis with the angular velocity u while the gear blank rotates about 

( 2 ) ( 2 ) 

the ^2 with the angular velocity w . Axes and intersect 

each other and form the angle 90® ^2 ~ ^ 2 * ^2 dedendum angle 

( 2 ) 

for gear 2. Axis X^^ is perpendicular to the generatrix of the root cone 
of gear 2. The coordinate system shown in Fig. 2.2 is rigidly connected 
to the housing of the gears and will be used for the analysis of conditions 
of meshing of the gears. 

Considering the generation of the pinion we use the following 

(F) 

coordinate systems: (i) which is rigidly connected to the generating 

surface (ii) which is rigidly connected to the frame of the 

xf lu 

cutting machine and (ili) which is rigidly connected to the pinion (gear 
(Fig. 2.3). Axes X^^^ and Z^ do not intersect but cross each other; AE^ 
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and ALj^ are the corrections of machine~tool settings which are used for the 

improvement of meshing of the gears. In the process of generation the 

generating surface rotates about the -axis with the angular velocity 

CP) 

oi while the gear 1 blank rotates about the -axis with the angular 
velocity Axes and Z^ form the angle 90° - , where A is 

the dedendum angle of gear 1; axis is perpendicular to the generatrix 

of the root cone of gear 1 . 


3. 6ENEKATING TOOL SURFACES 

The tool surface is a cone and is represented in the coordinate system 

as follows (Fig. 2.1) 

s 


x^j^ 

s 


r^^ ^cotij;^^^ - u.cosTjj^^^ 
c J c 



u.sinip^^^sinS. 
J c J 

s 

II 

u.sinU'^^ ^cos6, 
J c J 

1 


L' 


(j = F, P) 


(3.1) 


where u^ and 9^ are the surface coordinates. 

The coordinate system ^ (j = F, P) is an auxiliary coordinate 
system which is also rigidly connected to the tool (Fig, 2.1,b). To 
represent the generating surfaces and S in coordinate system we 


use the following matrix equation 
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c 


s 

y(j) 
^ c 

cs 

,(J) 

s 



z^J> 

c 


s 

_ 1 


1 

^ 1 


1 

0 

0 

0 








s 

0 

cosq . 
J 

+ slnq . 
J 

+b .sinq . 
3 J 





S 

0 

+sinqj 

cosq , 
J 

b . cosq , 
3 j 


z^j> 



s 

0 

0 

0 

1 


1 


Here: and are parameters which determine the location of the tool in 

the coordinate system \ Henceforth, the upper sign corresponds to the 
generation of a left-hand spiral bevel gear that is shown in Fig. 2.1 and 
lower sign for a right-hand spiral bevel gear. 

Equations (3.1) and (3.2) yield 


c J c 


c 

, 0 ) 

c 

,(J) 


1 , 0 ). 


u.sinif; "^sin( 6.+q . ) + b.sinq. 
J c 11 11 

u ^cos(0 ,+ q , ) + b.cosq. 
J c j j J J 


(3.3) 


where j = (F, P). 

The unit normal to the generating surface Z 


j 


(j = F, P) is represented 
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by 


N 


(j) 


n 


(j) _ 




3r«> Srp> 
, where — X r-- 

~c ae, 9u. 

J 3 

(j) 


Using Eqs. (3.3) and (3.4) (provided u^sinij^^-' 0), we obtain 


(3.4) 


nO) _ sini|j^^ ^i^^ ^ + cosij;^^ ^ [sin(0,+q . ) ^ + cos( 0, - q.)k^^^] 
" T„„ 3 3 ~c 33 ~c 


c ~c 


(3.5) 


4. EQUATIONS OF MESHING BY CUTTING 

Generation of We derive the equation of meshing of the generating 

surface 2 and the gear tooth surface 2 using the following procedure: 

F i 

Step 1: First, we derive the family of surfaces that we represent in 

the coordinate system Such a family is generated while the 

(F3 (1) 

coordinate system is rotated about the X' -axis (Fig. 2.3). We 

(F) 

recall that the generating surface Ep is rigidly connected to . The 

(F) ( 1 ) 

coordinate transformation in transition from S to S is represented by 

c m 

the following matrix equation 


x^^> 


'x(F)' 

m 


c 

m 

m c 





m 


c 

_ 1 


1 


(4.1) 


Here (Fig. 2.3): 
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1 


0 


0 


0 


,„(1)(F) 
m c 



0 


0 


cos(})p sin4)p 0 

-sin4>p cos(})p 0 

0 0 1 


where <})_ is the angle of rotation about the -axis. 

^ m 

Using Eqs. (4.1), (4.2) and (3.3), we obtain 



m 




m 

m 

= 

UpSin’J^^^^sinXp + bpSin(qp 

F 


m 


(f) 

UpSinij;^ cosTp + bpCos(qp 

? <1>F^ 


1 


1 



where 

p F + ^F 


(4.2) 


(4.3) 


The upper sign corresponds to the right-hand spiral bevel pinion. 

Equations (4.3), with parameter 4>p fixed, represent a single surface of the 
family of generating surfaces. 

Step 2 ; The unit normal to the generating surface Up may be represented in 
the coordinate systems as follows: 



where N 


( 1 ) _ 


3r 


( 1 ) 




v.m 


ae. 


3r(l> 

~m 


(4.4) 


9 


We may also use an alternative method for the derivation of the unit 
normal. This method is based on the matrix equation. 


[„»>) . 


Matrix may be determined by deleting the 4^^ column and rows in 

m c 

matrix (4.2). 

(F) 

The column matrix [n ] is given by vector equation (3.5). 

c 

After transformations, we obtain 


[.<«] . 

m 


(F) 

cosip^ 


,(F) 

C0S1|J COSTp 

c f 


Step 3: We derive the equations of the relative velocity, v , as 



follows: 


«-m -vm 


( F) ( 1 ) 

where ' is the velocity of a point N on surface and v is the 
~m ^ 

velocity of the same point N on surface 

(F) 

Vector V is represented by the equation 
~m 


v<^> . 

~m ~m 


Here (Fig. 2.3): 
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(4.9) 


[co<f>l - 

m ** 


- 0 ) 


(F) 


Vector is represented by equation (4.3), 

Equations (4.8) and (4.9) yield 




m 


.(1) 


m 


-y. 


( 1 ) 


m 


(4.10) 


Gear 1 rotates about the Z - axis with the angular velocity 

(Fig. 2. 3). Since does not pass through the origin of the 

~ m 

coordinate system we substitute by an equal vector which passes 

m ~ 

through and the vector moment represented by 


(1) 

0 0, X 03 ^ 

m h .s. 


Then, we represent v^^^ as follows 

-m 


m h ~m 


~ m ~m m 


(4.11) 


Here: 
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= 

m n 


LsinAj^ 




-al. 


, and 

m 


= ^co 


(1) 


sin(y^ - Aj) 


cos(Yj - Ap 


(4.12) 


where L = 0,M. 

h 

Equations (4.11) - (4.12) yield 


‘'’1' *1 - 


- 0 ) 


(1) 


(1) 


-m 


sln( Yj 

„(1) 




m 


.o) u(i) 


-m 


m 


-m 


cos(Yj - Aj) 
_( 1 ) 


m 


(1) 


0) 


(1) 




LsinA, 


j(i) ^(1) 

~m 


-AE, -al. 


sin(Yj - A^ 0 


cos(yj - Ap 




+AE^)cos(Yj - Aj) 

“ LsinApcos(Yj ~ A^) ~ (z^^^+ AL^ sin(Y^ 


Al) 


(y. 


( 1 ) 


m 


+ AEp sin(Y^ - A^) 


(4.13) 


( FI ) 

The final expression for ' is 

-m 
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+ AE,)cos(y, “ Ai ) 

m i i i 


m 


“ LsinApcos(Yj - A^) 


m 




- + AL^)sin(Yj ~ Aj)| 


+ AEpsin(Yj - A^) 


(4.14) 


Step 4; The equation of meshing by cutting is represented by 


n 


( 1 ) 


~in 


,<"> . 0 


(4.15) 


Using equations (4.15), (4.14), (4.6) and (4.3), 
-Up + [r^^^cotijj^^^ - LsinAj^ - AL^ tan(Yj - 


AEj^ [slnij/^^ - cosiJJ^^^tan(Y “ Aj^)cosTp] = 

fl(up, 9p, ^) = 0 (4.16) 

Here; 


,(F) 


(F). 


simp sin(qp ^ (j)^) + cosi|; sin0p 


“f 1 - sln(Yp - 


cos(Yj - 


1 - V _ 

J 


we obtain 


(F) ^ 

Aj)]cos,j;^ t 


_ W' ' 

“f1 (1) 

(A) ' 

This equation relates the generating surface coordinates (up, 0p) with the 
angle of rotation (4>p)« 
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Generation of ^ 2 * By using a similar procedure we may obtain the equation 
of meshing for surface and surface as follows 

Step 1 : Equations of the family of generating surfaces represented in 

( 2 ) 

the coordinate system S are: 

m 


ra 


r^^^cot;!;^^^ - u cosih^'^^ 
c ^c P ^c 

m 

= 

u sinilj^^^slnT + b sin(q x (b ) 
P ^c P P ^p + ^p' 

4 '’ 


u sin^^^^^cosT + b cos(q r 4) ) 
i_ P c p p '^p + ^p'^ J 


( 2 ) 

Step 2: The unit normal to ^ is represented in ^ 

p 


(4-17) 


by the column matrix 




c 

c 

cos4>^^^ 

c 


sinT 

P 


COST 

P 


where 


(4.18) 


T 

P 


6 + q + 4> 

P P P 


The upper sign corresponds to the left-hand spiral bevel gear and the lower 


to the right-hand spiral bevel gear. 

(P) (2) (P2) (2) 

Step 3: The velocities v , v and v are represented in ^ as 

-111 ~m -“m 

follows 


m 


14 


0 


m 


(p) 


= 0) 


.( 2 ) 


in 


,( 2 ) 


m 




= - 0 ) 


( 2 ) 


y^2^cos(Y2 - Ap 

+ LsinA9)cos(Y, 
in L £. 

y^^^sin(Y, - A ) 
m / z 


A2> 


sin(Y^ 
m ^ 




Step A; The equation of meshing 


„<» . v“’^> . 0 




yields the following equation 




[up - + LsinA2)cosip^^^]sinTp + 


(?) (?) ‘VZ " sin(Y2-A2> 

bpSinCqp 4>p)sini^^ + bpCOSij;^ sin0p cos(.y^-K^) 


= 0 


Here: 


m 




(P) 


P2 (2) 

03^ 


(4.19) 



(4.20) 

(4.21) 


(4.22) 
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5. ORIENTATION OF THE PINION CRADLE 


Henceforth, we will consider two auxiliary coordinate systems and 
S^, that are rigidly connected to the gear and pinion cutting machines, 
respectively (Fig. 5.1). 

Figure 2.2 shows coordinate system that is rigidly connected to the 

( 2 ) 

gear cutting machine and to coordinate system . Coordinate system 

is also rigidly connected to the housing of the gear train and the meshing 

of the generated pinion and gear will be also considered in coordinate 

system S^. Axis is the instantaneous axis of rotation of the pinion and 

the gear - the pitch line (the line of tangency of the pinion and gear 

pitch cones). The origin 0^ of coordinate system coincides with the 

( 2 ) 

points of intersection of 3 axes: (Fig. 2.2) and (Fig* 2.3). 

( 2 ) 

Here: X is the axis of rotation of the gear cradle; Z is the axis of 

m ° 2 

rotation of the gear being in mesh with the generating gear Zp and Z^ 
is the axis of rotation of the pinion being in mesh with the generating 
gear Zp and the gear member. 

Figure 2.3 shows coordinate system that is rigidly connected to the 

pinion cutting machine and coordinate system The origin 0, of 

m n 

coordinate system coincides with the origin 0^ but the orientation of 
system with respect to represents a parameter of the machine-tool 
settings (proposed by Litvin, 1968). This parameter is designated by e 
in Fig. 5.1. 

The coordinate transformation in transition from to is based on 
the following considerations. 

(i) Consider two auxiliary coordinate systems, S and S, , that are 

a D 

rigidly connected to systems and respectively (Fig. 5.1, a and 

Fig. 5.1,b), axes Z and z. coincide with the pinion axis. Initially 

a D 

coordinate system coincides with and system with S^. 
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(ii) Assume now that coordinate systems and are rotated about the 
Z -axis, the axis of the pinion, through the angle £ (Fig. 5.1,c). 

di 

This angle determines the orientation of coordinate system with 
respect to S^, or the orientation of the pinion cutting machine with 
respect to the gear cutting machine. Matrix [M^^] represents the 
coordinate transformation in transition from to and is represented by 
the following equation 


t^fhl = ^^fa’ f\bl ^h’ 

0 
0 
0 
1 

(5.1) 

6. PLANE OF NORMALS 

It will be proven below that the generating surfaces Zp and Zp contact 
each other at a point that moves in the same plane, II . The generated 
pinion and gear surfaces, and Z^, also contact each other at every 
instant at a point that coincides with the point of tangency of surfaces Z^ 
and Zp. However, we have to emphasize that surfaces Zp and Z^^ 
(respectively, Zp and Z^) are in line contact and the instantaneous line of 
contact moves over the contacting surfaces. The special property of the 
method developed for generation of spiral bevel gears is that the contact 
point moves in a plane (plane II) that is rigidly connected to the fixed 


2 2 
cosecos sin Yj^ 


sinecosYj^ cosYj^sinYj^(l ~ cose) 


-sinecosY 


1 


cose 


sines inY, 


2 2 

cosYjSinYj^ (1 ~ cose) -sinesinY|^ cosesin Y^ + cos Yj^ 
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coordinate system S^. Also, the common normal to the contacting surfaces 
lies in plane J[ and, as it will be shown later, performs a parallel motion 
in the process of meshing. Plane II is called the plane of normals and it 
is determined as the plane that passes through the instantaneous axis of 
rotation, Z^, and the normal N to the generating surface Zp (Fig. 6.1). 

The above normal passes through point N that is the main gear contact 
point. 

Figure 6.2 shows the orientation of coordinate system with respect 

to coordinate system S^. Origin 0^ coincides with origin 0^ and axis 

coincides with the Z - axis. The coordinate transformation in transition 

n 

from to is given by the matrix 



cosn 

-sinri 

0 


iJ-fnl - 

sinfi 

cosri 

0 

(6.1) 


0 

0 

1 



The determination of angle 0 is based on the following considerations: 
(1) unit vector i can be represented in coordinate system S- by the 

-v^n I 

following matrix equation 




cosn 

-Sinn 

0 


" 1 “ 

[if”']- lifnllin' 


Sinn 

cosn 

0 


0 



0 

0 



0 


cosri 






Sinn 






_0 
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( 6 . 2 ) 


(p) 

(2) The unit noraal to the generating surface is represented as 

follows 


, (P), ,, (2)„ (2) 


(6.3) 


Here (Fig. 2.2) 


[L 



COSA 2 0 -sinA2 

0 1 0 

sinA2 0 COSA 2 


(6.4) 


Unit vector n is given by column matrix (4.18). Equations (6.3), (6.4) 
•^xn 

and (4.18) yield 


(«f>i 


(P) (p) 

sin^l^^ cosA^ - cos4'^ cosTpSinA2 


(P) 

cos4>^ sinTp 


( P) (P) 

sinip '^sinA- + cosij^ cost^cosA, 

C Z c r Z 


(6.5) 


(3) 


Vectors if”^ and nf^^ are mutual perpendicular, i.e., il'^^ • n^ ^ - 0. 


~f 


This yields that 


tanr) 


(P) 

tanij^^ 'COSA 2 - cosTpSinA2 
sinT„ 


( 6 . 6 ) 


Here: 


6p+qj 


+ <Pr 


19 


with the upper sign for the left-hand gear, 

(P) 

ijj' = a„ for the convex gear tooth side, 

C if 

(P) 

^ 180® - ap for the concave gear tooth side 
where ap is the angle of the cutter blade. 

Equation (6.6) determines the angle of orientation n of the plane of 
normals (Fig. 6.2). 

7. PERFORMANCE OF PARALLEL MOTION OF A STRAIGHT LINE PROVIDED BY TWO 
RELATED ELLIPSES. 

An important part of the proposed approach is a new technique directed 
at the performance of a parallel motion of a straight line provided by two 
related ellipses. By using this technique it becomes possible to provide a 
parallel motion for the common normal to the gear-pinion tooth surfaces. 
This motion is performed in a plane (the plane of normals) that is rigidly 
connected to the gear housing and has the prescribed orientation. 

It is well known that a translational motion of a straight line may be 
performed by a parallelogram linkage (Fig. 7.1, a). Consider that a 
straight line slides by its points A and C along two circles of equal 
radii. Vectors v. and v^ which represent the velocities of points A and C 
of the moving straight line are equal. The moving straight line AC being 
initially installed parallel to the center distance OD will keep its 
original direction in the process of motion. 

The discussed principle of translational motion of a straight line may 
be extended for the case where the straight line slides along two mating 
ellipses (Fig. 7.1,b). These ellipses have the same dimensions and 
orientation and again the velocity vectors and y^ are equal. Consider 
that the moving straight line is initially installed parallel to OD where 0 
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and D are the foci of symmetry of the ellipses. Then, with the 

moving straight line will keep its original direction in the process of 
motion. 

Figure 7.1(a) and Fig. 7.1(b) show a translational motion of a segment 
of a straight line (AC) with constant length. A more general case is 
represented in Fig. 7.1(c). The straight line slides over two ellipses 
whose dimensions and orientation are different. The length of segment AC 
which slides along the ellipses is changed in the process of motion. The 
problem is how to provide a parallel motion of the the moving straight 
line. We call this motion a parallel one because the straight line has to 
keep its initial parallel to line OD where 0 and D are the foci of symmetry 
of two mating ellipses. Unlike the cases which are shown in Fig. 7.1(a) 
and Fig. 7.1(b) the motion of straight line is not translation because the 
velocities y^ and v^ of the tracing points A and C are not equal. The 
distance between the sliding points A and C is changed in the process of 
motion. It will be proven that the required parallel motion of straight 
line may be performed with certain relations between the dimensions and 
orientation parameters of two guiding ellipses which are shown in Fig. 

7.1(c). 

Consider that an ellipse (Fig. 7.2) is represented in coordinate 
system by the equations 

= 0, = BpCOsyp, = bp siniip (7.1) 

where ap and bp are the lengths of the semlmajor and seinlmlnor axes, respec- 
tively, D is the origin of coordinate and the symmetry focus of the ellipse 
parameter Up determines the location of a point on the ellipse. 
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Coordinate system 

has the 

same origin as S and the 
a 

orientation of 

with respect to is given by angle 

6p. To represent the ellipse in 

coordinate system we use the following matrix equation. 


n 

[M ][r^^^] 
na a 




(7.2) 

where 

“1 0 

0 

0 



na 

0 cos6p 

-sin6p 

0 


(7.3) 


0 sin6p 

cos6p 

0 



1 

0 

0 

1 

0 

1 




Matrix equation (7.2) yields 


(P) A (P) X U 4 X • 

x^ = 0, = apCos6pCosyp - bpSinfipSinpp 


(P) 

= apSin6pCosyp + bpCos6pSinyp 


(7.4) 


The S3rmmetry focus of the mating ellipse is point 0 given by 
coordinates: =0, ^ ^3 “ ^3 algebraic 

values). The orientation of the ellipse is given by angle 6^ (Fig. 7.2). 
Equations of the mating ellipse are represented in coordinate system as 
follows 


(F) (F) 

x^ = 0, y^ = apCos6pCosyp - bpSin6pSinyp + C 2 , 

(F) 

Zn = apSin6pCosyp + bpCOs6pSinyp + (7.5) 
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Consider point A of ellipse 1 is determined by parameter UpQ* Point C 
of ellipse 2 is determined with parameter Henceforth we will consider 

such a motion where Up - Up^ = Up - UpQ = U* Line AC is drawn through 
points A and C. Our goal is to provide that line AC will be parallel to 
the center distance DO for any value of the motion parameter U and perform 
a parallel motion. The above-mentioned goal can be achieved with certain 
relations between the ellipse parameters Sp, bp, 6p, Sp, bp and 6p. We 
start the derivation of these relations by considering the vector equation 


DA(U) + AC(U) = D0(U) + OC(U) 


(7.6) 


Vector equation AC(u) = ^DO is satisfied with any value of U if the 

parallel motion of line AC is provided. (Where X is the constant required 

to make the two vectors of equal magnitude.) 

Equation (7.6) yields 


(OC-DA)-j^ (X-l)(DO-j^) 
(OC-DA)*k (X-l)(DO*k) 


cos q 
sin q 


(7.7) 


Here: j and k are the unit vectors of axes Y and Z , q is the angle 

in ~n n n 

formed by axis Y^ and vector DO (Fig. 7.2). For the further 

transformation we will represent U-n Vi-r. ss follows 

P F 


Up = Upo + u, Up = Upo + u 


(7.8) 


Equations (7.7), (7.4), (7.5) and (7.8) yield 


(bj^j^sinq - aj^j^sinq - b2j^cosq + a2^cosq)cosU 


+ (-bj^ 2 ®^^^ b.,.,cosq + a, ,sinq - a,,cosq)sinU = 0 


"22 


12 “ 


’ 22 ' 


(7.9) 
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Here: 


^11 

^12 

®21 

^22 


apCosfipCOSUpQ - 
apCos6pSinypQ + 
apSin6pCosypQ + 
apSin6pSinypQ - 


bpSln6pSinypQ 

bpSln6pCOSypQ 

bpCos6pSinypQ 

bpCos6pCosypQ 


(7.10) 


Coefficients b^j, bj 2 » b 2 j and b 22 have similar expressions. 

Equations (7.9) must be satisfied for any value of the motion parameter y . 
This means that equation (7.9) will be satisfied for any value of y if the 
following two equations are satisfied simultaneously. 


bjjSinq - b2jCOsq - a^^sinq + a2j 
bj^2®^^9 ” b22COsq - aj^2®^^9 "*■ ^22 


cosq = 0 
cosq = 0 


(7.11) 

(7.12) 


Equations (7.10), (7.11) and (7.12) yield a system of two pseudo-linear 
equations in two unknowns (cosyp^ and sinyp^) : 


apSin(q-6p)cosypQ - bpCOs(q-5p)sinypQ - d^ 
bpCos(q-6p)cosypQ + apSin(q-6p)sinypQ = d2 


(7.13) 

(7.14) 


Here: 


dj = apSin(q-6p)cosypQ - bpCos(q-5p)sinypQ (7.15) 

d2 = apSin(q-6p)sinypQ + bpCOs(q-<Sp)cosypQ (7.16) 
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The solution of equations (7.13) and (7.14) for cosy^Q and siny^^ is 


as follows 


cosy^o 


apdiSinCq-Sp) +bpd 2 COs(q- 6 F) 

V 


apSin^(q-6F) + bFCOs^(q-6 ) 


siny 


FO 


3^(12 s in (q-Sp) -bFd^cos(q- 6 F) 
apSin^Cq-fip) +b^cos^(q-6^) 


F' ■ "F 

Equations (7.13) and (7.14) yield 


2.2. X N ^ v2 2. . . j2 ^ .2 

apSin (q-6p) + bpCOS (q-6p) = d^ + 


(7.17) 


(7.18) 


(7.19) 


Equations (7.19) relates 3 parameters of the mating ellipse: a^, b^, 

and (Parameters Sp, bp, 6p of the first ellipse and q are considered 

as given). Thus Eq. (7.19) can be satisfied with various combinations of 
ap, bp and 6p. 

Example; a^ = 2, b^ = 1.25, 6^ = 20% q = 270% = 3, y^^ = yp^ = 0. 

F 

Then we obtain: a_ = 2.2753, 6 = 34.3113® . These ellipses are shown in 

r r 

Fig. 7.2. 

A linkage which may perform the described parallel motion of line AC is 

shown in Fig. 7.3. Link 1 is in contact with two perpendicular slots and 

DN = a , DM = b where 2a and 2b are the lengths of the axes of ellipse 
P P P P 

which is traced out by point A. Similarly, link 1' is in contact with two 
other perpendicular slots, ON' = a„ and OM' = b , point C traces out 

r r 

ellipse 2. The orientation of a pair of slots with respect to the other 

one is determined by the angle (6„ ~ 5 ) (see Fig. 7.2). Links 1 and 1' 

j; p 

rotate with the same angular velocity w = Line AC will perform a 

parallel motion if the ellipse parameters satisfy equation (7.19). 
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8. GEAR MACHINE-TOOL SETTINGS 


Input Data * The following data are considered as given to set up the gear 
machine-tool settings . 

the teeth number of the gear 

Y 2 » the gear pitch angle (Fig. 2.2) 

the gear dedendum angle (Fig. 2.2) 

the mean spiral angle (Fig. 2.1) 

L , mean distance of the pitch cone (Fig. 2.2) 

W , the point width of the gear cutter (Fig. 8.1) 

the blade angle of the gear cutter (Fig. 8.1) 


Gear cutting ratio . The gear cutting ratio represents the ratio between the 
angular velocities of the cradle and the generated gear and is designated 
by 


“p2 


Ui 


(p) 


0 ) 


( 2 ) 


( 8 . 1 ) 


( 2 ) 

Henceforth, we will consider the two coordinate systems, and 

shown in Figure 2.2. Axis is the pitch line of the mating spiral bevel 
gears. It is also the instantaneous axis of rotation of the pinion and the 
gear that transform rotation with constant angular velocity ratio. Assuming 
that the generating cone surface, the gear tooth surface, 

pinion tooth surface, are in continuous tangency at every instant, it 
is required that the instantaneous axis of rotation by the gear cutting must 
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coincide with the pitch line, axis Z 


Thus 


f 


. ^(P) _ ^(2) . ,, 

~ ~ ~ ~ r 

This means that vectors and k, are collinear. 

~ ~ I 

Vectors and are represented in coordinate system 

~ ~ ~ X 

as follows 



0 


■-sin(Y 2 - ^2)“ 
0 


0 

1 _ 


cos(Y 2 - ti ^)_ 


sinA^ 


COSA2 


Equations (8.2) and (8.3) yield that 


+ o)^^^sin(Y2 - A^) 
(2l 

'^cos(Y2 - ^2) 


= tanA2 


Equation (8.4) results in that 

o ,( P > ^^"^2 

“p2 ^(2) C 0 SA 2 


( 8 . 2 ) 




(8.3) 


(8.4) 


(8.5) 
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Main Gear Contact Point . The main gear contact point is the center of the 


bearing contact on the gear tooth surface. The location of this point can 

be represented by the surface coordinates of the tool cone 0 and u . 

P P 

We can vary the value of u^ by keeping 0^ fixed and obtain the desired 

location of the main gear contact point. 

"k 

We start with the case when coincides with the pitch point 

k 

(Figure 8.2). In this case N on the cutter surface is the contact point 
of the pinion and gear tooth surfaces since the normal at this point 
intersects the instantaneous axis of rotation, Z^. The basic relations 
for the left-hand gear for this case are as follows: 


0 = 90° -B +q =T +q 
P P P P P 


b 

P 


LcosA2Sin(0p - q^) 
sin0 


LcosA_cosB 
^ P 

cos (Bp - qp) 


( 8 . 6 ) 


(8.7) 


d LcosA^sinq LcosA,.sinq 
J = 2 P = 2 > 

2 -* 


( 8 . 8 ) 


sin0 


cos (Bp - qp) 


a = -^ + LsinA«tanijj 
2 2 c 


(8.9) 


, , d LcosA„sinq 

r<P) = -£ - a = ^ 

’^c 2 ^ 


£ _ 


sin0 


^ - LsinA2tanifj^^^^ 


( 8 . 10 ) 


* 

u = 
P 


sini|; 


(P) 

- + asini|< 

(p) c 


( 8 . 11 ) 
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For right-hand gear we need the following changes: 


e* = 270® + 6 
P P 



(8.6a) 


LcosA„ cos3 LcosA- cos3 

V, = 2 p ^ 2 p 

p sin(9* - 180®) cos(8 - q ) 
p P p 


(8.7a) 


d LcosA„ sinq 
P ^ 2 

2 sin(6* - 180®) 
P 


LCOSA2 

cos(3p 



(8.8a) 


Let us now consider point N whose surface coordinates of the tool cone 
* 

are 0 and 
P 


(■o') fP') * * * 

u = 0^‘^'^N = 0^ •'n + N N = u + Au 

P P P 


( 8 . 12 ) 


Point N will become the point of contact if the normal vector of this point 
passes through M^. This requirement will be observed if the cradle is turned 
at a certain angle 4>p that is shown in Figure 8.3. We can determine angle (})p 
using the equation of meshing of the generating gear and the generated gear , 
equation (4.22), together with equations (8.5), (8.6), (8.7), (8.8), (8.9), 
(8.10), (8.11), and (8.12). This yields 


f2(AUp, (j>p) = 


Au cos(3 + <}> ) ± LcosA^sinil^ ^^^sinij) + LsinA-cos;)^ [cos3 “ cos (3 + <}>J] 
p p p Z c p Zc p PP 


(8.13) 
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The upper sign corresponds to the left-hand spiral bevel gear and the lower 
to the right-hand spiral bevel gear. Based on trigonometry transformations 


a closed form solution for (J)^ is obtained. That is 


*^p -A ± VA^ + 

2 


(8.14) 


A = cosA„sinifJ^^^ - sinA^cosiJ^ '■^‘^sinB + sinB 

^ c c o L 


(p)„.. 


B = — cos3 - sinA^cosij; ^ cosg 
L p 2 p 


C = sinA^cosij; 

2 l 


For gear concave side using the same equations except 


a = 2 


- LsinA2tanijj^ 


(8.15) 


, X d LcosA_sinq ,, . 

P 


(8.16) 


we can get the same result. 

Equation (8.14) determines the main gear contact points for both convex 


side and concave side. In general, AUp > 0 for gear convex side, and Au^ < 0 
for gear concave side. The absolute value of AUp for gear convex side is 
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greater than that of Au^ for gear concave side. 


Direction of the Contact Path . The direction of the contact path may be 
determined by the direction of the tangent to the contact path at the main con- 
tact point. This can be done by differentiating equation (4.22). In the pro- 
cess of motion the contact normal performs a parallel motion. Thus parameter 

T is a constant and this yields that 
P 

dT =0 d6„ = -d(t) (8.17) 

P P P 

Taking into account equations (8.17) and (8.5), we obtain after the differen- 
tiation of equation (4.22) the following equation 


du 


a A sinT + b^cos(q + 4> )sini() 
d())p p P '^p ^p' ^c 


(P) 


b cost); 
P c 


(P) 


COS0 tanA2 = 0 


The upper sign corresponds to the left-hand gear and the lower to the right-hand 

* * 

gear. At the main contact point we have that 6 = 6 , cos0 = sin(6 - a ). 

P P p p T * 

LCOSA2COSB 

b_ = -- .0 Then we obtain 

p cos(3p - qp) 


du Lcos3 P ... 


sinA2sin(3p - qp)cosij^^^^j = 0 (8.18) 


This yields 

E 

tan q ® - — 
^P D 


(8.19) 
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Here: 


D = 


du 

I 

d({) 


sinB slnx 
P P 


■t LcosA„sini|^^^^cosB sind) 

2 p 


LsinA-cosd'^^^cos^3 

2 p 


du 


,(P) 


(P), 


E = cos3 sinT + LcosA„sinij^ ‘^'^cos3 cosd) - LsinA«cosTjJ^^'^sin3 cos3 
dd)„ P P 2 P ^P 2 ^p 1 


There is a particular case when the direction of the contact path is 

du d 

parallel to the root cone and thus = 0. For this case, however, 

d(f)p 2 ^ 

becomes substantially larger than LC 0 SA 2 . More favorable ratio for 2 LcosA 

du 2 

can be obtained by decreasing the value of , but a more inclined path 

^p du 

of contact will occur. A reasonable value of is about -0.5 for the 

aq>p 

left-hand gear convex side, or -0.2 for the left-hand gear concave side. 


Because both sides of the gear are cut by duplex method, the values of 
du du 

-r^ for both sides are not independent, and we have to choose the value of -r^ 
^‘Pp ^p 

for both sides by a compromise. 


y. BASIC PRINCIPLES FOR GENERATION OF CONJUGATE GEAR TOOTH SURFACES 


IMTRODDCTION 

The method for generation of conjugate spiral bevel gear tooth surfaces 
is based on the following principles: 

(1) Four surfaces-two generating surfaces (Zp and Ep) and pinion and gear 
tooth surfaces (Ej^ and E 2 ) in tangency at every instant. The ratio of 
the angular velocities in motion of surfaces Xp> Ep» Ej ®nd ^2 ™nst satisfy 
the following requirements: (i) the above-mentioned surfaces must be in 

continuous tangency and (ii) the generated pinion and gear must transform 
rotation with zero kinematical errors, and (iii) Point N of contact of both 
pairs of contacting surfaces - Ep Ep, E^ E 2 ” moves in plane II (Fig. 9.1) 
in the process of meshing and the common normal to above-mentioned surfaces 
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performs parallel motion in plane H, 


Figure 9.1 shows the drawings that are represented in the plane of 

normals, II* Point A is the point of intersection of the gear head-cutter 

axis with plane II and A is the initial position of the normal to surfaces 

Ep and E^. Point D is the point of intersection of the gear cradle axis 

with plane H. Simultaneously point D is the point of intersection of 

pinion and gear axes. Axes of the gear cradle and gear head-cutter are 

parallel (the head-cutter axis is not tilted with respect to the axis of 

( 2 ) 

the cradle). The cradle axis, , is perpendicular to the gear root cone 

but it is inclined with respect to the plane of normals, II, 

Point C (Fig. 9.1) is the point of intersection of the pinion tool cone 

axis with plane II . Surface of the pinion tool cone will be in contact 

with surface (and ^ 2 ) common normal to surfaces Ep and E 2 passes 

through point C. Simultaneously, the pinion surface, Ej^ , will be in 

contact with E (and E and E«) at point N if the equation of meshing for 

surfaces and Ep is satisfied at N . 

Point 0 that is shown in Fig. 9.1 is the point of intersection of the 

pinion cradle axis with plane H. Point 0 does not coincide with point D of 

the gear cradle axis. We recall that specific machine tool-settings for 

the pinion, AE^ and AL^ , (Fig. 2.3) have to be used to provide conjugate 

pinion and gear tooth surfaces. This method for generation provides the 

collinearity of vectors DO and NA. We recall that the pinion cradle axis, 

is perpendicular to the pinion root cone but it is inclined with 
m 

( 2 ) 

respect to gear cradle axis, , and plane H. The orientation of axes 

X^^^ and X^^^ with respect to each other and plane II depends on gear and 
m m 

pinion dedendum angles, Aj^ and The pinion cradle axis and the pinion 

head-cutter axis are parallel (the pinion cutter is not tilted with respect 
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to the pinion cradle). 

Figure 9.1 shows the instantaneous positions of the contact normal N 
and points A and C. In the process of meshing the normal performs a 
parallel motion in plane II while point A (and respectively C) traces out an 
ellipse with the ellipse symmetry center D (symmetry center 0, 
respectively). The dimensions and orientation of the gear ellipse centered 
at D are known since the gear machine settings are given. The dimensions 
and orientation of the ellipse centered at 0 must be determined to provide 
the desired parallel motion of the contact normal. Then it becomes 
possible to determine the pinion machine-tool settings. 


10. SATISFACTION OF THE EQUATIONS OF MESHING 

Four surfaces “ ^2 ” contact within the 

neighborhood of the instantaneous contact point if they have a common 
normal and the following equations of meshing are satisfied at the point of 
contact. 

(1) The equation of meshing for surfaces and is represented as 
follows (Fig. 9.1): 


v<*2) . „ . , ,»)) . „ = . 0 


( 10 . 1 ) 


Here: oi 


( 12 ) ( 1 ) ( 2 ) ^ ( 1 ) ^ ( 2 ) ^ , ... 

I (i\ = 0 ) - 0 ) where o) and o) are the angular velocities 

( 12 ) 


of the pinion and the gear, respectively. Vector w represents the 

velocity in relative motion. Vector ^ is the position vector of contact 

( 12 ) 

point N and N is the contact normal. Vector w must be directed along 
the pitch line to provide transformation of rotation with the prescribed 
angular velocity ratio. 
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(2) The equation of meshing of surfaces 2p and ^2 represented by the 
equation 

= 0 ( 10 . 2 ) 

CP21 (Pi (2) 

Equation (10.2) is satisfied ±f dy .= Oi ' - (ly is collinear to the 
unit vector a of the pitch line. This requirement is observed already 
since the gear cutting ratio is determined with equation (8.5). 

(3) The equation of meshing of tool surfaces 2p and Sp is represented as 
follows : 




N 




X r 


(N) 


) + (DO X I • N = 0 


(10.3) 


(F) 

Deriving equation (10.3) we substitute the sliding vector ^ that 
passes through 0 (Fig. 9.1) by an equal vector that passes through D and 
the vector moment represented by 


m = DO X OJ 


(F) 


In the process of motion the contact normal keeps its original direction that 
is parallel to DO. Thus equation (10.3) yields 




[( 0 )' 


,(P) 


- J^>) 


.(N) 


N] = 0 


(10.4) 


Equation (10.4) may be interpreted kinematically as a requirement that 
vectors r^^^ and N must line in the same plane, the plane of normals 

n. This means that 




i - 0 

~n 


(10.5) 
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where is the unit vector of the -axis that is perpendicular to plan II. 

It was mentioned above that the proposed method for generation is based 
on the parallel motion of the contact normal that slides along two related 
ellipses. It was assumed that this motion is performed with the following 
conditions (see Section 7): 




dt dt dt 


( 10 . 6 ) 


where VL and 1-L are the ellipse parameters. 

Let us prove that the requirement can be observed if 

I = I and equation (10.5) is satisfied. 

(F) (P) (1) (2) 

Vectors ^ and ^ are directed along the - and x axes and 

^ -s. mm 

equation (10.5) may be represented as follows 


/f) . I - a,'"'* . 1 


(10.7) 


(F) (P) (F) (P) 

Since w' ' = w = oo we obtain that projections of vectors o) and oj 

on the normal to the plane are equal. 

This yields that 


^ IH 

dt dt 


Here: 


^ = o)(i;^^^ • i ) = • i ) 

at ~m '-n '-m ^n 


[i^^^ • 1 


( 10 . 8 ) 
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dy 

Vectors and — — represent the angular velocities of links 1 and 1 

dt dt 

7.3) in their rotational motions. Equations (10.8) also yields that 


/ 


(Fig. 


id) . i = id) . i (10.9) 

^ XII 


Requirement (10.9) can be satisfied with a specific orientation of the 
pinion cradle coordinate system, Sj^, with respect to gear cradle coordinate 
system, S^. This orientation may be achieved with a certain value of angle 
e (Fig. 5.1,c) that may be determined by using equation (10.9). The matrix 
representation of equation (10.10) is given as follows: 


11 0 



II 0 01 



( 10 . 10 ) 


Here.matrix Is given by Eq. (6.4); matrix is the transpose 

matrix of given by equation (6.1); matrix [L^^l is the sub-matrix of 

[Mj^l given by equation (i>.l); matrix i® represented as 

follows (Fig. 2.3); 


r cosA^ 0 




-sinAj^ 0 


sinAj^~ 

0 

cosAj 


Equation (10.10) yields 


( 10 . 11 ) 
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( 10 . 12 ) 


cos cose - sinetanri 


cosA^ - sinYj^sin(Y^ - A^ 
cos(Yj - A^) 


Equations (10.12) provides two solutions for angle £ and the smaller value 
of s is to-be chosen. 


11. EQUATIONS OF SURFACE TANGENCY AT THE MAIN CONTACT POINT 

We consider that the coordinates of the main contact point N on surface 

( P) 

(Fig. 8.2) and the direction of the surface unit normal n at point N 

are known. We have to provide that surfaces E , E and E will be in 

P Z 1 F 

tangency at the chosen point N. Surfaces E^ and E^ are already in tangency 
at N since the equation of meshing of E^ and E^ is satisifed at this point 
(see Section 8). Then we have to provide the tangency of surfaces and 
Zp, Ep and E^ and and 
Equations of Tangency of Ep and E p 

The tangency of generating surfaces 2p and at point N is provided if 
the following equations are satisfied at N: 






.(F) 


, AE^, AL^ = r 


= JP) 


( 11 . 1 ) 


.(P)/ 


„(F) 


n-^(0p. 4>p, qp) = n 


(P) 


( 11 . 2 ) 


(p) (p) 

Here r and n are the given position vector and surface Ep unit 
normal. Vector Eqs. (11.1) and (11.2) may be considered in any coordinate 
system, for instance, in system * Vector Eq. (11.1) provides three 

independent scalar equations but Eq. (11.2) provides only two ones since 
|n| = 1. Parameters Up and 6p are the surface Ep coordinates for the point 
of tangency. Parameter c()p, is the angle of cradle turn for the installment 
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of the generating surface Zj, and it may be chosen of any value, including 

(})p = 0. Parameters and determine the location and orientation of the 

(1) ( F) 

pinion head-cutter in coordinate system S ' (Flg» 2.1); lii^ ' is the blade 

m ^c 

angle (Fig. 2.1); AE^ and AL^ are the machine-tool settings that have been 

(F) 

shown in Fig. 2.3; r^ is the radius of the head-cutter circle obtained by 

the intersection of the head-cutter surface with plane = 0 (Fig. 2.1). 

m 

We may represent vector equations of tangency (11.1) and (11.2) as 
follows 




(11.3) 

(11.4) 


( 2 ) 

Matrix is represented by the following equation (Fig. 11.1) 




cosA^ 0 -sinA2 LsinA2CosA2" 


1 0 


sinA 2 0 C 0 SA 2 Lsin A 2 


0 0 


(11.5) 


Matrix is the transpose matrix of matrix [Mp^] that is represented by 

equation (5.1). 

Matrix [M ] is represented as follows (Fig. 2.3) 
m n 
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[M = 

m h 


cos 


0 -sinAj^ LsinAj' 


-AE. 


sinAj 0 cosAj -AL^ 


0 0 


( 11 . 6 ) 


(P) (P) 

The column matrix and [n ] are represented by equation (4.17) and 

(4.18). After matrix multiplications we obtain 


(1) (2)j ^ 

m m 


11 


^3 

^14 


21 

CM 

^23 

CM 

(11.7) 

'31 

^32 

'^33 

^34 


1 

0 

0 

1 



Here; 


hi = 

cose cos (Y j 

- A^)cos(Yj 

+ A2) + sin(Yj - Aj)sin(Yj + A2) 

^12 “ 

-sinecos(Y ^ 



h3 “ 

-cosecos(Yj 

- Aj)sin(Y 

j + A2) + sin(Yj - Aj)cos(Yj + A2) 

h4 = 

LsinA^ + LsinA2 bj^j 


^21 

sine cos(y ^ 

+ A 2 > 



cos e 
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= - sin£sin(Yj + 


\>2i^ - - AE^+ LsinA2 b2j^ 

b^j^ = -cosesin(Yj^ ■ Aj^)cos(Yj + + cos(Yj^ - Apsln(Y]^ + A 2 ) 

b ^2 “ sinesin(Yj - A^ 

b ^2 = cos£sin(Yj “ A^sin(Yj^ + + cos(y^ - A^cosCy^ + 

’"SA “ ^31 

Equations (11.3) and (11.4) are used for the determination of pinion 
machine-tool settings. 

12* PINION MACHINE-TOOL SETTINGS 

The pinion machine-tool settings are represented by the following 
(F) 

parameters: r^ - the radius of the head-cutter circle measured in plane 

X =0; b„ and q„ - parameters that determine the location of the 
m r r 

head-cutter in plane X = 0 with <|)_ = 0 (see Fig. 2.1); 

® ' (F) 

- parameter that determines the blade angle; m^j^ = ~7yT ~ the ratio by 

^ iiT ^ 

cutting; AE^^ and ALj^ - corrections of machine-tool settings that are shown 
in Fig. 2.3. 

The determination of the pinion machine-tool settings is based on the 
equations that relate the parameters of the gear and pinion ellipses, 
equations of tangency of the pinion and gear tooth surfaces at the main 
contact point and the equation of meshing by pinion cutting. These 
settings must be determined for the following 4 cases that are represented 
in the plane of normals, by Fig. 12.1 and Fig. 12.2, respectively. Figure 
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12.1(a) corresponds to the case where the gear is left-hand, the gear tooth 

(F) (P) 

contacting surface is convex, and CN > AN. Figure 12.1(b) 

corresponds to the case where the gear is left-hand, the gear tooth 

( F) (F) 

contacting surface is concave, r < r and CN < AN. 

c c 

Drawings of Fig. 12.2 correspond to the cases where the gear is 

right-hand and the gear tooth contacting surface is convex (Fig. 12 .2, a) 

(P) (F ) 

and concave (Fig. 12. 2, b). The difference between r and r with other 

c c 

parameters as given determines the required dimensions of the contacting 
ellipse. 

Determination of the Cutting Ratio 

It is mentioned above that the proposed method for generation provides: 
(i) the simultaneous tangency of four surfaces: Sp, and ^ 2 * where 

Zp and Zp are the gear and pinion generating surfaces and Z^^ and Z 2 are the 
generated pinion and gear tooth surfaces; (ii) the generating gears 
(cradles) are rotated with the same angular velocities i.e. co^ = by 
This yields the following equation for the pinion cutting ratio 


0) 03 m. 


Here: 


>2 equation 8.5) 


^'2 


12 (2) N, 

0 ) ' 1 


where N 2 and are the numbers of gear and pinion teethe 
The final expression for the pinion cutting ratio is 
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( 12 . 1 ) 


Nj^sinY2 
”Vl N2COSA2 


Determination of 

The procedure of computations of b„ is based on the following steps; 

(i) determination of semiaxes of the gear cradle ellipse, Sp and bp; (ii) 
determination of the parameter of orientation of the gear cradle ellipse, 

6p; (iii) determination of the parameter of orientation of the pinion 
cradle ellipse, 6p and, (iv) determination of b^. 

Step 1: Determination of semiaxes of the gear cradle ellipse. 

The gear cradle ellipse is obtained by intersection of the cylinder 
of radius bp (Fig. 2.1) with the plane of normals. The semiminor axis of the 
cradle ellipse is the cylinder radius bp. The semimajor axis is represented 
by the equation 


.( 2 ) 


( 12 . 2 ) 


i 

~m ~n 

( 2 ) 

where the unit vector i is collinear to the cylinder axis, i is the 

unit vector that is perpendicular to the plane of normals. The scalar 

( 2 ) 


product i 


m 


• i^ may be represented as follows 


i(2) 


i 

~n 


[1 0 0][L^(2)^][LjJ 



(12.3) 


Here, the column matrix represents the unit vector i^ in coordinate system 

( 2 ) ( 2 ) 

S and the row matrix represents the unit vector i in S \ 
n ~m 

Equation (12.3) yields 
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i 


( 2 ) 


= cos^cosri 


and 


(12.4) 


a 


P 



COSA2COSTI 


(12.5) 


Step 2 ; Determination of orientation of the gear cradle ellipse. 

o 

We designate with bp the unit vector of the minor ellipse axis and with 
6p the angle that is formed by axis and the major ellipse axis (Fig. 
12.3). Vectors b° and i^^^ are mutual perpendicular and this yields that 


ro -I 


[1 ooia«)^iiLjj 


-sin6- 


siniSpSinhcosA^ + cosSpSinA^ = 0 


cos 6, 


( 12 . 6 ) 


and 


tanA„ 

tan&, r— 

r stnri 


(12.7) 


Step 3 ; Determination of orientation of the pinion cradle ellipse. 

The unit vector of minor axis of the pinion cradle ellinse, b°, (Fig. 12.3) 

is perpendicular to the cylinder axis of radius b-. Thus 

r 



The matrix representation of equation (12.8) is as follows: 
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r 0 


11 0 0UL<1) I -sin(^ 


= 0 


cos 




(12.9) 


Equation (12 .9) yields 


tan6„ = 4^ (12.10) 

F B 

where 

A =» cosesinYj^ “ cosYj^tan(Yj “ Aj) (12.11) 

B = cosnsine + sinTi[ cose cosy sinY^tan(Yj “ A^)] (12.12) 

Step 4 ; Determination of 

It is easy to prove that the ratio between the axes of both ellipses is the 
same since 


.1 ,^. 1 < 2 ) . 1 

a„ -m -n -m .^n 

t P 

and 

• i = • i (see Eq. (10.9)) 

~m ~n <N^m •^n 

Thus : 


cosA^cosri ’ cos^cosn 


(12.13) 
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The parameters of both ellipses have been represented by equation (7.19) as 
follows 


apSin^Cq - 6^) + b^cos^Cq - 6^) = 

Equations (7.15), (7.16) and (12.12) yield that 


(12.14) 


d 


2 

1 




sin^(q-6 ) 

cos (q - 6p) + 2 2~ 

cos A 2 COS ri 


(12.15) 


where q is the angle that is formed between vector DO and axis (Fig. 
7.2). We recall that vector DO is collinear to the contact normal. Using 
equations (12.14) and (12.15), we obtain 




2 2 2 2 

cos (q-6p)cos A 2 COS o + sin (q - <Sp) 
_ _ _ _ 

cos (q - 6p)cos A 2 COS 0 + sin (q - 6^) 


1/2 


(12.16) 


Determination of q ^ 

The main idea of determination of qp is based on identification in 
plane = 0 of projections of two vectors: OC (Fig. 12.1 and Fig. 12.2) 

and a that is the unit vector of O^^^O^^^ (Fig. 12.4). We recall that 

'w ms 

points O' ' and 0 lie on the axis of the cradle and O' and C lie on the 

axis of the pinion head-cutter (Fig. 12.4). Points 0 and C are the points 

of intersection of pinion cradle axis, and the pinion head-cutter 

m 

axis, 0^ 'Cy with the plane of normals, !!• Plane IT has been shown in 

s 

Fig. 6.1. 


Unit vectors i , a and OC lie in the same plane. 
~m ~ 

unit vector of axis 

m 

Thus: 


Here i 


( 1 ) 


~m 


is the 
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(12.17) 


[a OC] = a 

~ ~TH ~ 


X OC) = 0 
~m 


The determination of the is based on equation (12.17) and the 
computational procedure may be presented as follows: 

Step 1 : Representation of vector a in coordinate system 

Figure 12.5(a) and Fig. 12.5(b) show the orientation of unit vector a 
in plane in two cases where the generating gear is left-hand 

(generated gear is right-hand) and generating gear right-hand (generated is 
left-hand), respectively. Vector a is represented by the following column 
matrix: 


0 


a 


sinqp 

cosqp 


(12.18) 


Angle qp is measured clockwise. This angle is negative if measured 
in opposite direction as shown in Fig. 2.1. 

The generating gear generates the member-gear with the same direction of the 
spiral (the member-gear is down with the respect to the generating gear as 
it is shown in Fig. 2.2). The generating gear generates the pinion with the 
opposite direction of the spiral (the pinion is up with the respect to the 
generating gear as it is shown in Fig. 2.3). 


Step 2 : Representation of vector OC in coordinate system S^. Using 

equations (7.5) and (12.13) we may represent vector OC as follows: 
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0 


OC = apCos6pCosypQ - bj,sin6pSinUpQ 


apSin6pCOsypQ + bpCosSpSlnyp^ 


ap cosSpCOsypQ - sinSpSinypQCOsncosA^ 


sin6pCosypg + cos6pSinypQCOsricosA2 


(12.19) 


The sub-script "FO” in indicates that the initial position of the 

contact normal (at the main contact point) is considered. 

Step 3: Representation of vector OC in coordinate system 

m 

The coordinate transformation in transition from S to is 

n m 

represented by the product of matrices. 




11 

^12 

^13 

'21 

^22 

^23 

31 

®32 

^33 


Then vector OC may be represented as follows 


^ 12^2 ^IS'^S 


OC = ^ 22^2 ^ 23^3 


^ 32^2 


Step 4 ; Representation of the cross-product i^ x OC. 
It is easy to verify that: 
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0 


i X OC = 




^ 22^2 ■*■ ® 23^^3 


(12.22) 


Step 5 ; Determination of 

Equations (12.17), (12.18) and (12.22) yield 

-sinqp(a^2^2 ^33*^3^ cosqp(a22b2 + ^23^3^ “ ^ 


Thus 


tanqp 


^22^2 ^23^3 

®32^^2 ®33^3 


(12.23) 


Here: 


a 22 “ -sinEcosYjSinn + cos£cosD (12.24) 

a22 = -sinesinYj^ (12.25) 

a^2 “ sinntcosecosYj^sin(Yj - Aj) - sinyj^cos(Yj “ A^] 

+ cosDsinesin(Yj^ ~ Aj^) (12.26) 

a^j = cosYjCOs(Yj^ “ Ap + cosesinYjSin(Yj ~ A^) (12.27) 

Expressions of 62 and have been represented by the column matrix (12.19) 

in terms of UpQ* Our purpose is to derive relations between PpQ* V>pQ aud 

y^^^ and Here: y„_ is the parameter of the gear cradle axis for 

n n PO 
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m ■) (N ) 

point A (Fig. 9.1); ' atid are coordinates of the main contact 

point. Equations (7.17) and (7.18) yield: 


tan y. 


tan(q - 6 ) - cosncosA 


tan(q - 6^) +-^ cosncosA^ 


( 12 . 28 ) 


Here (see equations (7.15) and (7.16)) 


d^ tan(q - 6p) tanUp^ + cosncosA^ 
d7 ~ tan(q - 6p) - cosncosA^tanyp^ 


(12.29 ) 


The determination of parameter Pp^ is based on the following considerations 
(Fig. 12.6) 


0 A = 0 N + nA 
n n 


(12.30) 


Equation (12.30) yields 


- yj^^^= -AN • cosq 


-An * sinq 

n n 


Thus : 


,(A) . (N) 


(12.31) 


Using equations (7.4) and (12.5) we get 


(A) _ ^ / slnipccsyp^ 

^n P \ cosycosA^ P^ '^^PO 


(12.32) 
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(A) _ ^ r°"^pC°s^PO . „ . 

\ costicosA„ P®^^^PO 


(12.33) 


and are coordinates of the main contact point N in coordinate 

n n ^ 

system S^,q determines the orientation of the contact normal in plane 

X = 0. 

n 

Using Eqs. (12.31) - (12.33) we may determine parameter UpQ* 

Determination of Relation Between AE ^ and AL j^ 

Parameters AEj^ and AL^ of the pinion machine-tool settings have been 

shown in Fig. 2.3. Our goal is to prove that the proposed method for 

generation relates AEj^ and ALj^ in a certain way. 

Consider the drawing of Fig. 12.7. Point D coincides with the origins 

0^ of coordinate system S^, 0^ of coordinate system (Fig. 5.1) and the 

origin of coordinate system S^. Point D is also shown in Fig. 9.1. Point 

0 is the point of intersection of the pinion cradle axis with the plane of 

normals, II. Vector W must be collinear to the contact normal (Fig. 9.1). 

Figure 12.7 shows a plane that is drawn through points 0. 

Since vectors i^^\ and OD lie in the same plane, we have that 

~m m 


DO i ] = 0 
m ~m 


(12.34) 


Vector is represented by the column matrix (see Eq. 4.12)) 


O^^^D = = 

m m h 


LsinA, 


ae. 


al. 


(12.35) 


Vector DO is collinear to the unit contact normal n^ and is given by (see 
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Eq. 4.6)) 


[n - 
m 


slnip 


(F) 


(F) 

cosip sinT 
C r 


I (F) 

COS'p^ COSTp 


(12.36) 


Equations (12.34) - (12.36) yield 


LsinA 


1 


AE, 


AL, 


sintl)^ cosii> '^sinT- cosil; 'cost. 
c F V I 


(12.37) 


Equation (12.36) results in that 



tanTp 


where 


(12.38) 


"f ■ * "f *F 

(F) 

Determination of 4^. , t„ and AE, 

C r i 

(F) 

Parameter determines the blade angle of the pinion head-cutter. 

Parameter with the known value of q^ determines the cone surface 
parameter 0^. We may determine and by using the equations of 

tangency at the main contact point represented as follows: 
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(12.40) 


^(P) ^ /F) 

(P) _ (F) 

n =* n 


(12.41) 


Equations (12.40) and (12.41) provide the equality of position-vectors 
and unit surface normals, respectively at the main contact point for the 
generating surfaces Ep and Zp. Vector equations (12.40) and (12.41) 
provide three and two scalar equations, respectively, since | n'^ - |n^ '\ 

= 1. Using vector Eq* (12,41) we may determine parameters and T„, To 

C r 

determine AEj^ we need only one scalar equation from the three ones provided 
by vector equations but the remaining two scalar equations should be 
checked to ensure that they are satisfied with the determined machine-tool 
settings • 


13. INSTALLMENT OF MACHINE-TOOL SETTINGS 

The installment of the eccentric angle, EA, and cradle angle, CA, on 
the Gleason cutting machine No. 26 and No. 116 provides the required values 
of q^ and b^ (j = P, F), 

Consider that a left-hand gear is generated. Figure 13.1 shows two 
positions of the eccentric and its center: before and after the 

installment of the cradle angle, CA. The axis of the cradle is 
perpendicular to the plane of drawings and pointed to the observer. The 
axis of the head-cutter passes through the center of the cradle since the 
eccentric angle is not installed yet. Figure 13.2 shows two positions of 
the axis of the head-cutter: before and after the installment of the 

eccentric angle. It is evident that angles CA and EA are related by the 
equation 
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(13a) 


q. = CA + ^ - 90" (j = P. F) 

J 2 

= 00 sin ^ (13.2) 

2 c I 

where 00^ = k is the given constant value. Figures 13.3 and 13.4 
correspond to the installment of a right-hand gear (the generating gear is 
left-hand). Figure 13.3 shows the installment of the cradle angle (CA). 
Figure 13.4 shows the Installment of the eccentric angle and the relations 
between angular parameters CA, EA and . Equations (13.1) and (13.2) 
also for the case of generation of a right-hand gear but the value of q. is 
larger than in the case of generation of a left-hand gear (see Fig. 13.2). 

The installment of corrections of machine-tool settings Ae^ and AL^ are 
only applied for the generation of the pinion. These corrections are 
installed on the Gleason’s cutting machine by ’’the change of machine center 
to back” and ”the sliding base”. 

The correction Ae^ represents the shortest distance between the crossed 
axes of the pinion and the gear; AE^ is called ”the offset” in Gleason 
terminology (Fig. 13.5) . 

The correction Al^ is installed at the machine as a vector-sum of: (i) 

the change of machine center to back (CB) and (ii) the sliding base (SB) 
(Fig. 13.6a,b). The change of machine center to back is directed parallel 
to the pinion axis and the sliding base is directed parallel to the cradle 
axis. It is evident that 



CONCLUSION 


The authors proposed a method for generation of spiral bevel gears with 
conjugated gear tooth surfaces with the following properties: 

(i) The transformation of rotation is performed with zero or almost 
zero kinematical errors. 

(ii) The contact point of gear tooth surfaces (the center of instantaneous 
contact ellipse) moves in a plane (ir) of a constant orientation 
that passes through the pitch line of the gears. 

(iii) The normal to the contacting surfaces moves in the process of motion 
in plane it keeping its original orientation. 

(iv) It is expected that due to the motion of the contact ellipse along 
but not across the gear tooth surfaces the contact ratio will be 
increased and the conditions of lubrication improved. 

The authors developed a method of parallel motion of a straight line by 
two ellipses with related dimensions and orientation. 

A TCA program for the simulation of meshing and bearing contact for the 
proposed gears has been developed. The advantage of the proposed gearing is 
that the gears can be manufactured by the Gleason's equipment. 
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NUMERICAL EXAMPLE 


Basic Machine Settings 



Gear (LH) 

Pinion (RH) 

Both Sides 

Concave 

Convex 

# of Teeth 

N 

41 

10 

10 

Kean Spiral Angle 


35® 

— 

— 

Dedendum Angle 

A 

3° 53 ' 

1®41' 

1®41' 

Mean Cone Dist. 

B 

3,226 

— 

— 

Cutter Radius 

B 

3.957128371 

3.895686894 

4.034935029 

Cutter Width 

W 

0.08 

— 

— 

Blade Angle 

a 

20® 

16.7979304880® 

23.188994098® 

Machine Root Angle 


72.4097056667® 

12.0236276667® 

12.0236276667® 

Radial 

B 

3.37751754380 

3.31529203349 

3.45277453402 



73.6837239464® 

75.1337128849® 

71.9315208492® 

Machine Offest 

AE 

0 

0.005082532914 

-0.009816930736 

Mach. Ctr. to Back 

B 

0 

-0.00041258199 

0.00062147078 

Sliding Base 

Xs 

0 

-0.00008594705 

0.00012946173 

Ratio of Roll 

in 

0.973756061793 

0.237501478486 

0.23750147848 

Orientation - Angle 

£ 

! 

1 

1 

-0.26906022362® 

0.22644169356® 

Tilt 

B 

0 

0 

0 

Swivel 

j 

0 

0 

0 
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c. 
c. 
c. 
c. 
c. 
c. 
c. 
c 

C PURPOSE 
C 

c 
c 
c 

C NOTE 
C 
C 
C 
C 
C 
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* 
Vf 

MACHINE SETTINGS OF CONJUGATE SPIRAL BEVEL GEAR - 


AUTHORS: FAYDOR LITVIN 
WEI-JIUNG TSUNG 
HONG- T AO LEE 


vt 

Vc 


Vf Vf V? Vc Vf Vf Vc * Vf V? Vc Vr V? V? Vr * V? Vr Vc Vc Vr V? Vr Vf V? * ;V Vc ?V Vc * Vr >V Vr Vc V? Vc Vc Vc * V? V: V? Vc Vf Vt Vc Vc Vc Vc ?V Vr Vr V? Vr Vc V? Vc Vc Vc V? Vc 


To find the machine setting parameters for the generation of 
spiral bevel gears. 


This program is written in FORTRAN 77. It can be compiled by 
V compiler in IBM mainframe or FORTRAN compiler in VAX system. 

DESCRIPTION OF INPUT PARAMETERS 


c. . . 

JJ 


JJ=1 for left-hand gear, JJ=2 for right-hand gear 

c. . . 

TNI 

m 

teeth number of pinion 

c. . . 

TN2 

♦ 

teeth number of gear 

c.. . 

DIDG, DIMIN 

• 

dedendum angle of pinion (degree and arc minute, 

c. . . 



respectively) 

c. . . 

D2DG, D2MIN 

• 

♦ 

dedendum angle of gear (degree and arc minute, 

c... 



respectively) 

c.. . 

GAMADG 

• 

shaft angle (degree) 

c. . . 

BPDG 

• 

♦ 

mean spiral angle (degree) 

c. . . 

RL 

• 

mean cone distance 

c... 

W 


point width for gear finishing 

c.. . 

ALFAP 


blade angle of gear cutter (degree) 

c. . . 

RCFl 

• 

• 

pinion cutter radius for pinion concave side 

c. . . 

RCF2 

• 

# 

pinion cutter radius for pinion convex side 

c.. . 

DUPl 

• 

• 

a chosen value for pinion concave side to locate the 

c 



main contact point at the desired location 

c. . . 

DUP2 

• 

a chosen value for pinion convex side to locate the 

c 



main contact point at the desired location 

c... 

DUPFEP 

• 

a chosen value for pinion concave side to get the 

c 

r 



desired direction of the contact path 

C DESCRIPTION OF 

r 

OUTPUT PARAMETERS 

c. . . 

PHPDG 


blade angle of gear cutter (degree) 

c. . . 

PHFDG 


blade angle of pinion cutter (degree) 

c. . . 

RCP 


radius of gear cutter (measured on cradle plane) 

c... 

RCF 


radius of pinion cutter (measured on cradle plane) 

c.. . 

TPDG 


the difference between the sum of gear cone surface 

c 



angle coordinate and generating surface rotation 

c 



angle, and machine setting angle (degree) [see 

c 



equation (4.18) in report] 
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C. . . TFDG 

C 

C 

C. . . 

C... BP 
C... BF 

C... QPDG 
C. . . QFDG 
C... MP2 
C 

C. . . MFl 
C 

C... DEI 
C... DLl 
C 

C... FEEODG 
C 

C . . . DELTADG 

C 

C 

IMPLICIT 


: the difference between the sum of pinion cone 
surface angle coordinate and generating surface 
rotation angle, and machine setting angle (degree) 
[see equation (4.3) in report] 

: radial of gear cutter 
; radial of pinion cutter 
: machine setting angle of gear 
: machine setting angle of pinion 

; ratio of cradle angular velocity for cutting the 
gear and gear angular velocity 
: ratio of cradle angular velocity for cutting the 
pinion and pinion angular velocity 
; machine offset for cutting the pinion 
; vector-sum of (1) the change of machine center to 
back and (2) the sliding base 
: generating surface rotation angle at initial main 
contact point 

: rotation angle of frame h relative to frame f about 
Z1 

REAL*8 (A-H,0-Z) 


INPUT THE DESIGN DATA 


DATA JJ/1/ 

DATA TN1,TN2/10.D00,41.D00/ 

DATA D1DG,D1MIN/1.DOO,41.DOO/ 

DATA D2DG,D2MIN/3.DOO,53.DOO/ 

DATA GAMADG,BPDG/90.D00,35.D00/ 

DATA RL/3.226D00/ 

DATA W/0.O8DO0/ 

DATA ALFAP/20.D00/ 

DATA RCFl/3. 867079296 16D00/ 

DATA RCF2/4.05714282715D00/ 

DATA DUP1/0.042787628358701D00/ 

DATA DUP2/-0. 0377-63737441601673 IDOO/ 
DATA DUPFEP/-0. 526309620098410257D00/ 


CONVERT DEGREES TO RADIANS 

CNST=4 . D00*DATAN ( 1 . DOO) / 180 . DOO 
Dl= (DlDG+DlMIN/60, DOO) *CNST 
D2= (D2DG+D2MIN/60 . DOO) *CNST 
BP=BPDG*CNST 


CALCULATE GAMAl AND GAMA2 

GAMA=GAMADG*CNST 

RM12=TN2/TN1 

RM21=TN1/TN2 

GMA1=DATAN (DSIN (GAMA) / (RM12+DCOS (GAMA) ) ) 
GMA2=DATAN (DSIN (GAMA) / (RM21+DC0S (GAMA) ) ) 
IF(GMA1 .LT. 0.)THEN 
GMAl'180. *CNST+GMA1 
END IF 
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IF(GMA2 .LT. 0.)THEN 
GMA2=180.*CNST+GMA2 
END IF 

SUBSTITUTE SIN AND COS FUNCTION BY A SHORT NAME 

SNBP=DSIN(BP) 

CSBP=DCOS (BP) 

SND1=DSIN(D1) 

CSD1=DC0S(D1) 

SNR1D1=DSIN (GMAl-Dl) 

CSR1D1=DC0S (GMAl-Dl) 

SNR1D2=DSIN (GMA1+D2) 

CSR1D2=DC0S (GMA1+D2) 

SNR2D2=DSIN (GMA2-D2) 

CSR2D2*DC0S (GMA2-D2) 

CSD2=DC0S (D2) 

SND2=DSIN(D2) 

CSGM1=DC0S(GMA1) 

SNGM1=DSIN(GMA1) 

CSGM2=DC0S (GMA2) 

SNGM2=DSIN(GMA2) 

CALCULATE CUTTING RATIOS 

RMP2=SNGM2/CSD2 
RMF1=TN1*SNGM2/ (TN2*CSD2) 

11=1: THE PINION CONCAVE PART ANALYSIS, 11=2: THE PINION CONVEX PART 
ANALYSIS 

DO 10000 11=1,2 
IF (II .EQ. DTHEN 
RCF=RCF1 
DUP=DUP1 
ELSE 

RCF=RCF2 
DUP=DUP2 
END IF 

CALCULATE PHP 

IF (II .EQ .DTHEN 
PHP=ALFAP*CNST 
ELSE 

PHP= (180 . DOO-ALFAP) *CNST 
END IF 

SNPHP=DSIN(PHP) 

CSPHP=DCOS (PHP) 

CALCULATE FEEP AT INITIAL MAIN CONTACT POINT 

A=CSD2*SNPHP-SND2*CSPHP*SNBP+DUP*SNBP/RL 
IF(JJ .EQ. 2) THEN 
A=-A 
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END IF 

B=DUP*CSBP/RL-SND2*CSPHP*CSBP 

C=SND2*CSPHP*CSBP 

D=DSQRT (A*A+B*B-C*C) 

E=C-B 

F1=DATAN((-A+D) /E) 

F2=DATAN((-A-D)/E) 

IF (DABS (FI) .LT. DABS (F2) ) THEN 
FEE=2 . D00*F1 
ELSE 

FEE=2.D00*F2 
END IF 

CSFEE=DCOS (FEE) 

SNFEE=DSIN(FEE) 

CALCULATE TP 

IF(JJ .EQ. DTHEN 
TP=90 . D00*CNST-BP+FEE 
ELSE 

TP=270.D00*CNST+BP+FEE 
END IF 

SNTP=DSIN(TP) 

CSTP=DC0S (TP) 

CALCULATE QP 

IF (II .EQ. I) THEN 
IF(JJ .EQ. DTHEN 

D=DUPFEP*SNBP*SNTP+RL*CSD2*SNPHP*CSBP*SNFEE+RL*SND2*CSPHP 
// *CSBP*CSBP 

ELSE 

D=DUPFEP*SNBP*SNTP-RL*CSD2*SNPHP*CSBP*SNFEE+RL*SND2*CSPHP 
// *CSBP*CSBP 

END IF 

E=DUPFEP*CSBP*SNTP+RL*CSD2*SNPHP*CSBP*CSFEE-RL*SND2*CSPHP 
// *SNBP*CSBP 
QP=DATAN(-E/D) 

IF(QP .LT. O.DOO)THEN 
QP=I80.D00*CNST+QP 
END IF 
END IF 

IF(JJ .EQ. DTHEN 
SNQPFE'DSIN (QP-FEE) 

CSQPFE=DCOS (QP-FEE) 

THP=90 . D00*CNST+QP-BP 
ELSE 

SNQPFE=DS IN (QP+FEE) 

CSQPFE=DCOS (QP+FEE) 

THP=270 . D00*CNST-QP+BP 
END IF 

SNTHP=DSIN(THP) 

CALCULATE RCP 
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CM0=RL*CSD2*DSIN (QP) /DCOS (BP-QP) 

IF(II .EQ. DTHEN 

RCP=CM0-W/2 . 0D00-RL*SND2*SNPHP/CSPHP 
ELSE 

RCP=CM0+W/2 . D00-RL*SND2*SNPHP/CSPHP 
END IF 

CALCULATE RBP 

RBP=RL*CSD2*CSBP/DC0S (BP-QP) 


CALCULATE UP 

UP= (RCP*CSPHP/SNPHP+RL*SND2) *CSPHP+RBP* (SNPHP*SNQPFE-CSPHP* 
# SNTHP*SND2/CSD2) /SNTP 

CALCULATE RN 

AN=-SNPHP*CSD2+CSPHP*CSTP*SND2 

BN=CSPHP*SNTP 

RN=DATAN (AN/BN) 

SNRN=DSIN(RN) 

CSRN=DCOS (RN) 

CALCULATE DLTA 


A=CSGM1 

B=-SNRN/CSRN 

C= (SNGM1*SNR1D1-CSD2) /CSRIDI 
X1=-B+DSQRT (A**2+B**2-C**2) 

X2=-B-DSQRT (A**2+B**2-C**2) 

Y=C-A 

AE=X1/Y 

BE=X2/Y 

DLT2A=DATAN (AE) 

dlt2b=datan (be) 

I F (DAB S (DLT2A) . LT . DABS (DLT2B) ) THEN 
DLTA=DLT2A*2.D00 
ELSE 

DLTA=DLT2B*2.DOO 
END IF 

SNDLTA=DSIN (DLTA) 

CSDLTA=DCOS (DLTA) 

CALCULATE DETAP 

A01=SND2/CSD2 

B01=-SNRN 

DETAP=DATAN (aOI/BOI) 

CALCULATE DTAF 

A211=CSGM1*SNR1D1/CSR1D1-CSDLTA*SNGM1 

B21 1=- (SNDLTA*CSRN+SNRN* (CSDLTA*CSGMl+SNGMl*SNRlDl/CSRlDl) ) 
DETAF=DATAN (a21 1/B2I 1) 
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CSDTAF=DCOS (DETAF) 

SNDTAF=DSIN(DETAF) 

C 

C CALCULATE Q 
C 

C calculate normal vector at main contact point in F coordinate system 
C 

RNXF2=SNPHP*CSD2-CSPHP*CSTP*SND2 
RNYF2=CSPHP*SNTP 

RNZF2=SNPHP*SND2+CSPHP*CSTP*CSD2 

calculate normal vector at main contact point in N coordinate system 

RNYN2=-SNRN*RNXF2+CSRN*RNYF2 
RNZN2=RNZF2 

Q=DATAN2 (-RNZN2 , -RNYN2) 

QAP=Q-DETAP 
QAF=Q-DETAF 
SNQAP=DSIN(QAP) 

CSQAP=DCOS (QAP) 

SNQAF=DSIN(QAF) 

CSQAF=DCOS (QAF) 

CALCULATE BF 

BF=DSQRT ( (CSQAP**2*CSD2**2*CSRN**2+SNQAP**2) 

# / (CSQAF**2*CSD2**2*CSRN**2+SNQAF**2) ) *RBP 

CALCULATE PO 

RAP=RBP/ (CSD2*CSRN) 

calculate coordinate of main contact point in M coordinate system 

RMXM2=RCP*CSPHP/SNPHP-UP*CSPHP 
RMYM2=UP*SNPHP*SNTP-RBP*SNQPFE 
RMZM2=UP*SNPHP*CSTP+RBP*CSQPFE 

calculate coordinate of main contact point in F coordinate system 

RMXF2=CSD2*RMXM2-SND2*RMZM2+RL*SND2*CSD2 
RMYF2=RMYM2 

RMZF2=SND2*RMXM2+CSD2*RMZM2+RL*SND2*SND2 

calculate coordinate of main contact point in N coordinate system 

RMXN2=CSRN*RMXF2+SNRN*RMYF2 
RMYN2=-SNRN*RMXF2+CSRN*RMYF2 
RMZN2=RMZF2 
C 

AG=-RAP*SNQAP 

BG=RBP*CSQAP 

CG=RMYN2*DSIN (Q) -RMZN2*DC0S (Q) 

G1 1=-BG+DSQRT (AG**2+BG**2-CG**2) 
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G22=-BG-DSQRT(AG**2+BG**2-CG**2) 

G33=CG-AG 

P01=DATAN(G11/G33) 

P02=DATAN (G22/G33) 

IF (DABS (POD .LT. DABS (P02)) THEN 
P0=P02*2.D00 
ELSE 

P0=P01*2,D00 
END IF 

SNPO=DSIN(PO) 

CSPO=DCOS (PO) 

C 

C CALCULATE FO 
C 

AF=BF/ (CSD2*CSRN) 

RNF=BF**2 

RMF=AF**2 

DDll=RAP*SNQAP*CSPO-RBP*CSQAP*SNPO 
DD22=RAP*SNQAP*SNP0+RBP*CSQAP*CSP0 
SSS= (AF*SNQAF) **2+ (BF*CSQAF) **2 
CSFO= (AF*DD1 1*SNQAF+BF*DD22*CSQAF) /SSS 
SNFO= (AF*DD22*SNQAF-BF*DD1 1*CSQAF) /SSS 
F0=DATAN2 (SNFO , CSFO) 

C 

C CALCULATE QF 
C 

AA22=-SNDLTA*CSGM1*SNRN+CSDLTA*CSRN 

AA23=-SNDLTA*SNGM1 

AA32=SNRN*(CSDLTA*CSGM1*SNR1D1-SNGM1*CSR1D1)+CSRN*SNDLTA*SNR1D1 

AA33=CSGM1*CSR1D1+CSDLTA*SNGM1*SNR1D1 

BB22=AF*(CSDTAF*CSF0-SNDTAF*SNF0*CSRN*CSD2) 

BB33“AF*(SNDTAF*CSF0+CSDTAF*SNF0*CSRN*CSD2) 

UUU=- (AA22*BB22+AA23*BB33) 

DDD=AA32*BB22+AA33*BB33 
QFFE=DATAN2 (UUU , DDD) 

QF=QFFE 

SNQF=DSIN(QF) 

CSQF=DCOS (QF) 

SNQFFE=DSIN(QFFE) 

CSQFFE“DCOS (QFFE) 

C 

C calculate normal vector at main contact point in H coordinate system 
C 

RNXH2= (CSDLTA*CSGM1*CSGM1+SNGM1*SNGM1)*RNXF2-SNDLTA*CSGM1*RNYF2 
# +CSGMl*SNGMl* (1 .DOO-CSDLTa) *RNZF2 

RNYH2=SNDLTA*CSGM1*RNXF2+CSDLTA*RNYF2-SNDLTA*SNGM1*RNZF2 
RNZH2=CSGM1*SNGM1* (1 .DOO-CSDLTA) *RNXF2+SNDLTA*SNGM1*RNYF2+ 

// (CSDLTA*SNGMl*SNGMl+CSGMl*CSGMl) *RNZF2 

C 

C CALCULATE PHF 
C 

PHF=DAS IN (CSDl *RNXH2-SND1 *RNZH2) 

IF (II .EQ. 2) THEN 
PHF= 180 . DOO*CNST-PHF 
END IF 
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SNPHF=DSIN(PHF) 

CSPHF=DCOS(PHF) 

CALCULATE TF 

IF (II .EQ. 2) THEN 

TF=DATAN2 (-RNYH2 , -SND1*RNXH2-CSD1*RNZH2) 

ELSE 

TF=DATAN2 (RNYH2 , SND1*RNXH2+CSD1*RNZH2) 

END IF 

SNTF=DSIN(TF) 

CSTF=DCOS(TF) 

CALCULATE THF 

THF=TF+QFFE 

SNTHF=DSIN(THF) 

CSTHF'DCOS (THF) 

CALCULATE UF 

calculute coordinate of main contact point in H coordinate system 

RMXH2=(CSDLTA*CSGM1''CSGM1 + SNGM1*SNGM1)*RMXF2-SNDLTA*CSGM1>''RMYF2 

# +CSGM1*SNGMI* (1 .DOO-CSDLTA) *RMZF2 
RMYH2=SNDLTA*CSGM1*RMXF2+CSDLTA*RMYF2-SNDLTA*SNGMI''RMZF2 
RMZH2=CSGM1*SNGM1* (I .DOO-CSDLTA) *RMXF2+SNDLTA*SNGM1*RMYF2+ 

# (CSDLTA*SNGM1*SNGMI+CSGM1*CSGMI) *RMZF2 

UF= (RCF*CSPHF/SNPHF-RMXH2*CSD1+RMZH2*SND1-RL*SND1) /CSPHF 
CALCULATE DEl.DLI 

DE1=RMYH2-UF*SNPHF*SNTF+BF*SNQFFE 

DL1=DE1*CSTF/SNTF 

CONVERT RADIANS TO DEGREES 

PHPDG=PHP/CNST 

TPDG=TP/CNST 

QPDG=QP/CNST 

RNDG=RN/CNST 

DLTADG=DLTA/CNST 

DTAPDG=DETAP/CNST 

DTAFDG=DETAF/CNST 

QDG=Q/CNST i 

P0DG=P0/CNST 

F0DG=F0/CNST 

QFDG=QF/CNST 

PHFDG=PHF/CNST j 

TFDG=TF/CNST ' 

FEEDG=FEE/CNST 

WRITE (6, 60000) I 

WRITE(6, 60001) ' 

WRITE (6, 60002) I 

I 

I 

I 
I 
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IF(II.EQ.I) THEN 
WRITE (6, 60003) 

ELSE 

WRITE (6, 60004) 

END IF 

WRITE (6,60002) 

WRITE (6, 60001) 

WRITE (6,60000) 

WRITE (6, 90001) 

WR I TE (6 , 90002) PHPDG , PHFDG , RCP , RCF , TPDG , TFDG 

WR I TE (6 , 90003) RBP , BF , QPDG , QFDG , RMP2 , RMF 1 , DE 1 , DL 1 , FEEDG , DLTADG 

WRITE(6,1)RNDG 

WRITE (6, 20001) 

WR I TE (6 , 20002) DTAPDG , PODG , RBP , RAP 
WRITE (6, 20003) 

1 0000 WR I TE (6 , 20004) DTAFDG , FODG , BF , AF 


,'BP =' ,G18.12,15X, 'BF =',G18.12,/ 

,'QPDG =' ,G18.12,15X, 'QFDG =’,G18.12,/ 

,'MP2 =’ ,G18.12,15X, 'MFl =’,G18.12,/ 

,'DEl =' ,G18.12,15X, 'DLl =',G18.12,/ 

, 'FEEODG=' ,G18.12,15X, 'DELTA =' ,G18.12/) 

, 'ROTATION ANGLE OF RN =',G18.12/) 

, ' GEAR CRADLE ELLIPSE PARAMETER ' /) 

,'GEAR ELLIPSE ORIENTATION'' ,G24. 17,/ 

,' INITIAL POSITION'' ,G24. 17,/ 

, 'MINOR AXIS'' ,G24. 17, 5X, 'MAJOR AXIS' G24. 17/) 

, ' PINION CRADLE ELLIPSE PARAMETER ' /) 

, 'PINION ELLIPSE ORIENTATION' ' ,G24. 17, / 

,' INITIAL POSITION'' ,G24. 17,/ 

, 'MINOR AXIS'' ,G24. 17, 5X, 'MAJOR AXIS' ' ,G24. 17/) 

60000 FORMAT (IHl,' ') 

60001 F ORMA T ( 1 H ’ ^ * ) 

» Vf 

'* CONCAVE PINION PART ANALYSIS 

'* CONVEX PINION PART ANALYSIS 

'GEAR PARAMETER' ,27X, 'PINION PARAMETER') 

'PHPDG '' ,G18. 12, 15X, 'PHFDG '',G18.12,/ 

'RCP '' ,G18.12,15X, 'RCF =',G18.12,/ 

IH ,'TPDG '' ,G18.12,15X, 'TFDG '',G18.12) 


90003 FORMAT (IH 

# IH 

# IH 

# IH 

# IH 
1 FORMAT (IH 

20001 FORMAT (IH 

20002 FORMAT (IH 

# IH 

# IH 

20003 FORMAT (IH 

20004 FORMAT (IH 

# IH 

# IH 


60002 FORMAT (IH 

60003 FORMAT (IH 

60004 FORMAT (IH 

90001 FORMAT (IH 

90002 FORMAT (IH 

# IH 

# 

END 


,V.) 
* 1 ) 
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Q ^ ^ ^ ***«5V*»*****5C5C5V*5V**5t5C5V5V*5V**5V**?V*5V5V****yc5'cVc***i 

C . . . * 

C... * TOOTH CONTACT ANALYSIS OF SPIRAL BEVEL GEAR 

C... * 

C... * AUTHORS: FAYDOR LITVIN 

C... * WEI-JIUNG TSUNG 

C... * HONG-TAO LEE 

c . . . * 

C • • • 

Q ^ ^ ^ 'f * * * * * * ‘V * *'c Vf Vf >v Vf ?v y? Vc * y? y? yc y? y? y? y? yc y? >v y? y? y? yc yc y? y? yc y? yc y? y? y? : 

c 

C PURPOSE 
C 

C Tooth contact analysis of spiral bevel gears 


y? yr y? y? yc yc yc yc yc yc yc yc yc yc >v >v yc yc yc yc yc yc yc yc yc yc yc yc yc yc >v yc yc yc yc yc yc yc >v yc yc yc yc yc yc yc yc yc yc ?'c yc yc yc yc yc yc yc >v yc >v yc yc 


This program is written in FORTRAN 77. It can be compiled by 
V compiler in IBM mainframe. 

This program calls ZSPOW, a subroutine of IMSL package. 
DESCRIPTION OF INPUT PARAMETERS 


TNI 

TN2 

DIDG, DIMIN 

D2DG, D2MIN 

GAMADG 

RL 

PHPDG 

PHFDG 

RCP 

RCF 

TPDG 


JJ=1 for left-hand gear, JJ=2 for right*“hand gear 
11=1 for pinion concave side , 11=2 for pinion convex 
side 

teeth number of pinion 
teeth number of gear 

dedendum angle of pinion (degree and arc minute, 
respectively) 

dedendum angle of gear (degree and arc minute, 

respectively) 

shaft angle (degree) 

mean cone distance 

blade angle of gear cutter (degree) 
blade angle of pinion cutter (degree) 
radius of gear cutter (measured on cradle plane) 
radius of pinion cutter (measured on cradle plane) 
the difference between the sum of gear cone surface 
angle coordinate and generating surface rotation 
angle, and machine setting angle (degree) [see 
equation (4.18) in report] 

the difference between the sum of pinion cone 
surface angle coordinate and generating surface 
rotation angle, and machine setting angle (degree) 


c... 


[see equation (4.3) in report] 


c... 

BP 

: radial of gear cutter 


c... 

BF 

: radial of pinion cutter 


c... 

QPDG 

: machine setting angle of gear 


c... 

QFDG 

: machine setting angle of pinion [ALWAYS 

consider it 

c 

c 


as a POSITIVE input even it is negative 
the output of maching setting program 

gotten from 

i 

c... 

c 

MP2 

: ratio of cradle angular velocity for cutting the J 

gear and gear angular velocity 



non 


c. . . 

c 

MFl 


ratio of cradle angular velocity for cutting the 
pinion and pinion angular velocity 

c. . . 

DEI 

; 

machine offset for cutting the pinion 

c. . . 

c 

DLl 

• 

♦ 

vector-sum of (1) the change of machine center to 
back and (2) the sliding base 

c. . . 

c 

FEEODG 

• 

• 

generating surface rotation angle at initial main 
contact point 

c... 

c 

DELTADG 

• 

• 

rotation angle of frame h relative to frame f about 
Z1 

c. . . 

DFEPDG 

• 

# 

increment of FEPDG 


c 

IMPLICIT REAL*8(A-H,0-Z) 

INTEGER NS IG , lER , I TMAX , M , MLOOP , N 1 , N2 , N5 , NS IGl 
REAL*8 PAR5(5) ,X(5) , FNORM5 , WK5(75) ,MP2,MF1 
EXTERNAL FCNl,FCN2 
COMMON/Al/CNST 

COMMON/A2/DE1 , DLl , RCP , RL , RCF , PHFDG 

COMMON/A3/SNPHF , CSPHF , SNRIDI , CSR IDl , SNBP , CSBP , SNR 1D2 , CSR 1D2 
COMMON/aA/W , ECENDG , DI STEC 

COMMON/A5/SNGM1 ,CSGM1 ,GMAl , SNGM2,CSGM2,SND2,CSD2, SNDl ,CSD1 
C0MM0N/A6/SNTP , CSTP , TP , TPDG , SNTP , CSTF , TF , TFDG 
COMMON/A7/SNDLTA, CSDLTA, SNPHP ,CSPHP , SNR2D2 , CSR2D2 , DLTADG 
COMMON/ A 1 2/ UF , B F , RBP , UP , RMF 1 , RMP2 
COMMON/Al A/SNRN , CSRN , RNDG , RN 

C0MM0N/A18/SNQP , CSQP , QPDG , SNQF, CSQF , QFDG , QP , QF , QPO , QFO 
C0MM0N/A27/XF1 , YFl , ZFl ,XF2, YF2,ZF2 

COMMON/A32/AYN , AZN , B YN , BZN , CYN, CZN , DYN , DZN , RMYN , RMZN 
COMMON/A33/PHP 1 , PHF 1 , PHP2 , PHF2 , PHP IDG , PHF IDG , PHP2DG , PHF2DG 
C0MM0N/A34/FE1 , SNFEl , CSFEl , FE2, SNFE2 , CSFE2 , FEP , FE221 , FEl 1 1 , 

# FEPDG,FEFDG,F111DG,F221DG 

C0MM0N/A35/SNTHP , CSTHP , SNTHF , CSTHF , THPDG , THFDG 
COMMON /A36 /RKl 2 , RK22 , RKl 1 , RK2 1,UTX11,UTX12,UTY11,UTY12,UTZ11, 

// UTZ 1 2 , UTX2 1 , UTX22 , UTY2 1 , UT Y22 , UTZ2 1 , UTZ22 

COMMON/ A37/RNXF2 , RNYF2 , RNZF2 , RNXF 1 , RN YF 1 , RNZF 1 
COMMON/ A38 /DEF , S I GMDG , ALPHDG , AXI S A , AX ISB , ERROR 
COMMON/aAO/GMENRA, PMENR l , PMENR2 
COMMON/AA9/SNQAP , CSQAP , SNQAF , CSQAF 

COMMON/ A5 5 /DETAP , DETAF , SNDTAP , CSDTAP , SNDTAF , CSDTAF , DTAPDG , DTAFDG 
COMMON/A77/Q , QDG , PODG , FODG , PMIN , PMAX , FMIN , FMAX 

INPUT DATA 

JJ=1 
11=1 

TN1=10.D00 
TN2=A1.D00 
D1DG=1.D00 
D1MIN=A1.D00 
D2DG=3.D00 
D2MIN=53.D00 
GAMADG=90.DOO 
RL=3.226D00 
PHPDG=20.0D00 

PHFDG= 16 . 797930A880D00 ' 

RCP=3 . 83760775903D00 
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RCF=3.86707929616D00 
TPDG=53 . 0019620800D00 
TFDG=51.5519731A16D00 
BP=3.3775175A380D00 
BF=3.315292033A9D00 
QPDG=73 . 6837239A6AD00 
QFDG=77. 13175080A9D00 
MP2=.973756061793D00 
MF1=. 237501 A78A86D00 
DE1=0. 50825329 1A08D-02 
DL1=0. A03530720A3AD-02 
FEP0DG=-1. 9980379 199559D00 
DLTADG=- . 269060223623D00 
DFEPDG=1.00D00 
C 

C END OF INPUT DATA 
C 

DEF=0.00025D00 

N5=5 

NSIG1=12 

ITMAX=200 

RBP=BP 

RMP2=MP2 

RMF1=MF1 

C. . . 

C... CONVERT DEGREE TO RADIANS 

C. . . 

CNST=A . D00*DATAN ( 1 . DOO) / 180 . DOO 

GAMA=GAMADG*CNST 

QP=QPDG*CNST 

PHP=PHPDG*CNST 

PHF=PHFDG*CNST 

TP=TPDG*CNST 

TF=TFDG*CNST 

DLTA=DLTADG*CNST 

QF=QFDG*CNST 

FEPO=FEPODG*CNST 

FEFO=O.DOO 

Dl= (DlDG+DlMIN/60 . DOO) *CNST 
D2= (D2DG+D2MIN/60.D00) *CNST 
RM12=TN2/TN1 
RM21=TN1/TN2 

GMA1=DATAN (DSIN (GAMA) / (RM12+DCOS (GAMA) ) ) 
GMA2=DATAN (DSIN (GAMA) / (RM2 1+DCOS (GAMA) ) ) 

M= (360 . DOO/ (2 . 0D00*TN2) ) /DFEPDG 
CALL HEAD1(JJ,II) 

C. . . 

C... SUBSTITUTE SIN AND COS FUNCTION BY A SHORT NAME 

C. . . 

SND1=DSIN(D1) 

CSD1=DC0S (Dl) 

SNR1D1=DSIN(GMA1-D1) 

CSR1D1=DC0S (GMAl-Dl) 

SNR1D2=DSIN (GMA1+D2) 

CSR1D2=DC0S (GMA1+D2) 
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SNR2D2=DSIN (GMA2-D2) 

CSR2D2=DCOS (GMA2-D2) 

CSD2=DC0S (D2) 

SND2=DSIN(D2) 

CSGM1=DC0S(GMA1) 

SNGM1=DSIN(GMA1) 

CSGM2=DC0S (GMA2) 

SNGM2=DSIN(GMA2) 

SNPHP=DSIN(PHP) 

CSPHP=DCOS (PHP) 

SNQP=DSIN(QP) 

CSQP=DCOS(QP) 

SNTP=DSIN(TP) 

CSTP=DCOS(TP) 

SNDLTA=DSIN(DLTA) 

CSDLTA=DCOS (DLTA) 

SNQF=DSIN(QF) 

CSQF=DCOS (QF) 

SNPHF=DSIN(PHF) 

CSPHF=DCOS (PHF) 

SNTF=DSIN(TF) 

CSTF=DCOS(TF) 

TOOTH CONTACT ANALYSIS 

... DETERMINATION OF THE KINEMATICAL ERROR, BEARING CONTACT 

TP0=TP 

TF0=TF 

QP0=QP 

QF0=QF 

RINITP=M*DFEPDG 
FEP=-RINITP*CNST+FEPO 
C... INITIAL GUASS FOR ZSPOW 
X(1)=TF0 
X(2)=TP0 

X (3) =-RINITP*CNST+FEFO 

X(4)=0.0D00*CNST 

X(5)=0.0D00*CNST 

MLOOP=M*2+l 

MLOOPO=M+1 

DO 99 J=1,MLOOPO 

IF(JJ .EQ. DTHEN 

CALL ZSPOW (FCNl , NS IGl , N5 , ITMAX, PARS , X , FNORM5 , WK5 , lER) 
ELSE 

CALL ZSPOW (FCN2 , NS IGl , N5 , ITMAX, PARS , X , FNORMS , WKS , lER) 
END IF 

CALL ANGLE (X(D) 

CALL ANGLE (X (2)) 

CALL ANGLE (X (3)) 

CALL ANGLE (X (4)) 

CALL ANGLE (X(S)) 

FEF=X(3) 

FE2= (FEP-FEPO) /RMP2 
FE1=(FEF-FEF0) /RMFl 
FE111=X(4) 
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FE221=X(5) 

FE2PR=FE2-FE221 
FE1PR=FE1-FE111 
FEPDG=FEP/CNST 
FE1PRD=FE1PR/CNST 
FE2PRD=FE2PR/CNST 
FEP=FEP+DFEPDG*CNST 
99 CONTINUE 

FE1PR0=FE1PR 

FE2PR0=FE2PR 

RINITP=M*DFEPDG 

FEP=-RINITP*CNST+FEPO 

QP0=QP0 

QF0=QF0 

C... INITIAL GUASS FOR ZSPOW 
X(l)=TFO 
X(2)=TPO 

X (3) =-RINITP*CNST+FEFO 

X(A)=O.ODOO*CNST 

X(5)=O.ODOO*CNST 

MLOOP=M*2+l 

DO 88 J=l,MLOOP 

IF(JJ .EQ. DTHEN 

CALL ZSPOW (FCNl , NSIGl ,N5 , ITMAX, PARS ,X, FNORM5 , WK5 , lER) 

ELSE 

CALL ZSPOW (FCN2 , NS IGl , NS , ITMAX , PARS , X , FNORMS , WKS , lER) 

END IF 

CALL ANGLE (X(D) 

CALL ANGLE (X (2)) 

CALL ANGLE (X (3)) 

CALL ANGLE (X (4)) 

CALL ANGLE (X(S)) 

FEF=X(3) 

FE2= (FEP-FEPO) /RMP2 
FE1=(FEF-FEF0) /RMFl 
SNFE1=DSIN(FE1) 

CSFE1=DC0S(FE1) 

SNFE2=DSIN(FE2) 

CSFE2=DCOS (FE2) 

sntf=dsin(x(D) 

CSTF=DCOS (X(D) 

SNTP=DSIN(X(2)) 

CSTP=DCOS (X (2) ) 

FE111=X(4) 

FE221=X (S) 

FE2PR=FE2-FE221 

FE1PR=FE1-FE111 

ERROR' (FE2PR*3600 . D00-FE2PR0*3600 . DOO- (FE1PR*3600 . DOO-FEIPRO* 
# 3600.D00)*TN1/TN2)/CNST 

TFDG=X(1)/CNST 
TPDG=X(2)/CNST 
FEFDG=X(3)/CNST 
FEPDG=FEP/CNST 
F111DG=FE111/CNST 
F221DG=FE221/CNST 
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FE1PRD=FE1PR/CNST 

FE2PRD=FE2PR/CNST 

. . . DETERMINATION OF THE BEARING CONTACT 

• • • 

CALL FNORM (DN IXTH , DN IXFE , DN 1 YTH , DN 1 YFE , DN 1 ZTH , DN IZFE) 

CALL PNORM (DN2XTH , DN2XFE , DN2 YTH , DN2 YFE , DN2ZTH , DN2ZFE) 

CALL DUP(UPP,DUPTHP,DUPFEP) 

CALL DUF(UFF,DUFTHF,DUFFEF) 

CALL DR2 (UPP , DUPTHP , DUPFEP , DX2THP , DX2FEP , DY2THP , DY2FEP , 

# DZ2THP,DZ2FEP) 

CALL DR I (UFF , DUFTHF , DU FFEF , DXl THF , DX 1 FEF , DY 1 THF , DY I FEF , 

# DZITHF.DZIFEF) 

CALL PC (DN2XTH , DN2XFE , DX2THP ,DX2FEP , DN2 YTH , DN2 YFE , DY2THP , 

# DY2FEP,X12,X22,RK12,RK22) 

CALL PC(DN1XTH,DNIXFE,DX1THF,DX1FEF,DN1YTH,DN1YFE,DY1THF, 

# DY1FEF,X11,X21,RK11,RK21) 

CALL PD (DX2THP , DX2FEP , DY2THP ,DY2FEP , DZ2THP , DZ2FEP , X12 , 

// PDX12,PDY12,PDZ12) 

CALL PD (DX2THP , DX2FEP , DY2THP ,DY2FEP , DZ2THP , DZ2FEP , X22 , 

// PDX22,PDY22,PDZ22) 

CALL PD(DX1THF,DX1FEF,DY1THF,DY1FEF,DZ1THF,DZ1FEF,X11, 

# PDX1I,PDY11,PDZ11) 

CALL PD(DX1THF,DX1FEF,DY1THF,DY1FEF,DZ1THF,DZ1FEF,X21, 

# PDX21,PDY21,PDZ21) 

CALL UNIT(PDX11,PDY11,PDZ11,UNIT11,UX11,UY11,UZ11) 

CALL UNIT (PDX21 , PDY21 , PDZ21 , UNIT21 ,UX2I , UY21 , UZ21) 

CALL UNIT(PDX12,PDY12,PDZ12,UNIT12,UX12,UYI2,UZ12) 

CALL UNIT (PDX22 , PDY22 , PDZ22 , UNIT22 , UX22 , UY22 , UZ22) 

CALL UTRANl (UXl 1 , UY 1 1 , UZl 1 , UTXl 1 , UTYl 1 , UTZl 1) 

CALL UTRAN 1 (UX2 1 , UY2 1 , UZ2 1 , UTX2 1 , UTY2 I , UTZ2 1) 

CALL UTRAN2 (UXI2,UY12,UZ12,UTX12,UTY12, UTZ I 2) 

CALL UTRAN2 (UX22 , UY22 , UZ22 , UTX22 , UTY22 , UTZ22) 

CALL AX I S ( S I GMDG , AX I SA , AX I SB , ALPH , ALPHDG) 

CALL WRITE2 
FEP=FEP+DFEPDG*CNST 
88 CONTINUE 
STOP 
END 

• • • 

***** SUBROUTINE ANGLE ***** 

SUBROUTINE ANGLE (X) 

IMPLICIT REAL*8(A-H,0-Z) 

COMMON/ Al/CNST 

CNST2=2 . D00*CNST* 180 . DOO 

M-X/CNST2 

RM«M 

X=X-RM*CNST2 

RETURN 

END 

• • • 

^ ^ ^ ***** SUBROUTINE FCNl ***** 

# • • 

SUBROUTINE FCNl (X,F,N5, PARS) 
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IMPLICIT REAL*8(A-H,0-Z) 

INTEGER N5 

REAL*8 X (N5) , F (N5) , PARS (N5) 

COMMON/ Al/CNST 

COMMON/A2/DE1 , DLl , RCP , RL , RCF , PHFDG 

COMMON/ A3/SNPHF , CSPHF , SNR IDI , CSR IDI , SNBP , CSBP , SNR 1D2 , CSR 1D2 
COMMON/A5/SNGM1 , CSGMl , GMAl , SNGM2 , CSGM2 , SND2 , CSD2 , SNDl , CSDl 
COMMON/ A6/SNTP , CSTP , TP , TPDG , SNTF , CSTF , TF , TFDG 
COMMON/A7/SNDLTA, CSDLTA, SNPHP , CSPHP , SNR2D2 , CSR2D2 , DLTADG 
COMMON/A12/UF , BF , RBP , UP , RMFl , RMP2 

COMMON/Al 8/SNQP , CSQP , QPDG , SNQF , CSQF , QFDG , QP , QF , QPO , QFO 
COMMON/A27/XFl,YFl,ZFl,XF2,YF2,ZF2 

COMMON/A34/FE1 , SNFE 1 , CSFEl , FEZ , SNFE2 , CSFE2 , FEP , FE22 1 , FE 1 1 I , 

// FEPDG,FEFDG,F111DG,F221DG 

COMMON/A35/SNTHP , CSTHP , SNTHF , CSTHF , THPDG , THFDG 
COMMON/A37/RNXF2 , RNYF2 , RNZF2 , RNXF 1 , RNYF 1 , RNZF 1 
QF=QFO-X(3) 

SNQF=DSIN(QF) 

CSQF=DCOS (QF) 

QP=QPO-FEP 

SNQP=DSIN(QP) 

CSQP=DCOS (QP) 

FE111=X(4) 

FE22I=X(5) 

SN221=DSIN(FE221) 

CS221=DCOS(FE221) 

SN111=DSIN(FE111) 

CS111=DC0S(FE111) 

SNTF=DSIN (X(D) 

CSTF=DCOS (X(D) 

SNTP=DSIN(X(2)) 

CSTP=DCOS (X (2) ) 

THF=X(1)+QF 
THP=X (2) +QP 
THFDG=THF/CNST 
THPDG=THP/CNST 
SNTHF=DSIN(THF) 

CSTHF=DCOS (THF) 

SNTHP=DSIN(THP) 

CSTHP=DCOS (THP) 

UF1= ( (RCF*CSPHF/SNPHF-RL*SND1-DL1*SNR1D1/CSR1D1) *CSPHF*SNTF+ 

# BF* (SNPHF*SNQF+CSPHF*SNTHF* (RMFl-SNRlDl) /CSRlDl) -DEI* ( 

# SNPHF-CSPHF*CSTF*SNR1D1/CSR1D1) ) /SNTF 

XH1= (CSGM1*CSR1D1*CS1 1 1+SNGM1*SNR1D1) * (RCF*CSPHF/SNPHF-UF1* 

# CSPHF-RL*SND1) +BF* (CSGMl* (SNQF*SN1 1 l-CSQF*CSl 1 1* 

# SNR1D1)+SNGM1*CSQF*CSR1D1)+DL1*(SNGM1*CSR1D1-CSGM1* 

# SNR IDl *CS 111) -SNl 1 1 *DE1 *CSGM1-UF1 *SNPHF* (CSGMl* (SNTF* 

# SN111+CSTF*CS111*SNR1D1)-SNGM1*CSTF*CSR1D1) 

YH1=CSR1D1*SN1 1 1* (RCF*CSPHF/SNPHF-UF1*CSPHF-RL*SND1) -BF* ( 

# SNQF*CS111+CSQF*SN111*SNR1D1)-DL1*SNR1D1*SN111+DE1* 

# CSl 1 1+UF1*SNPHF* (SNTF*CS1 1 1-CSTF*SN1 1 1*SNR1D1) 

ZH1A= (CSGM1*SNR1D1-SNGM1*CSR1D1*CS111)*(RCF*CSPHF/SNPHF-UF1* 

# CSPHF) +RL*SND1*(SNGM1*CSR1D1*CS111-CSGM1*SNR1D1)-BF*(SNGM1 

# * (SNQF*SN1 1 1-CSQF*CS 1 1 1*SNR IDl) -CSGMl*CSQF*CSRlDl) + 

# DL1*(SNGM1*SNR1D1*CS111+CSGM1*CSR1D1)+DE1*SN111*SNGM1 
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ZH1=ZH1A+UF1*SNPHF* (SNGMl* (SNTF*SN11 1+CSTF*CS1 1 1*SNR1D1) +CSGM1* 
#. CSTF*CSR1D1) 

XF1= (CSDLTA*CSGM1**2+SNGM1**2) *XH1+SNDLTA*CSGM1*YH1+ (1 . ODOO- 
// CSDLTA) *CSGM1 *SNGM1 *ZH 1 

Y F 1 =-SNDLTA*CSGM 1 *XH 1 +CSDLTA* YH 1 + SNDLTA*SNGM1 *ZH 1 
ZF1= (1 . ODOO-CSDLTA) *CSGMl*SNGMl*XHl-SNDLTA*SNGMl*YHl+ 

# (CSDLTA*SNGM1**2+CSGM1**2) *ZHl 

XF2A=RBP* (CSPHP*SNTHP*SND2/CSD2-SNPHP*SNQP) * (CSPHP*CSR2D2 

# *CSGM2*CS221/SNTP+CSPHP*SNR2D2*SNGM2/SNTP-SNR2D2*SNPHP* 

# CSGM2*CS221*CSTP/SNTP+SNPHP*CSR2D2*SNGM2*CSTP/SNTP-SNPHP* 

# CSGM2*SN221) +CSPHP*CSGM2* (RCP*CSPHP+RL*SNPHP*SND2) * (SNTP* 

# SN221+CSTP*CS221*SNR2D2) 

XF2=XF2A+RBP*CSGM2* (SNQP*SN221+CSQP*CS221*SNR2D2) +RL*SND2* 
if (CSGM2*CSR2D2*CS221+SNGM2*SNR2D2)-RBP*CSR2D2*CSQP*SNGM2 

if +CSGM2*CSPHP*CSR2D2*CS221* (RCP*SNPHP-RL*CSPHP*SND2) -RCP* 

if CSPHP*SNGM2* (CSPHP*CSR2D2*CSTP-SNPHP*SNR2D2) -RL*CSPHP*SND2 

if *SNGM2* (SNPHP*CSR2D2*CSTP+CSPHP*SNR2D2) 

YF2=CSPHP*CSR2D2*SN221* (RL*CSPHP*SND2-RCP*SNPHP) +RBP* (SNPHP* 
if SNQP-CSPHP*SNTHP*SND2/CSD2) * (SNPHP*CS221-SNPHP*SNR2D2* 

if SN221*CSTP/SNTP+CSPHP*CSR2D2*SN221/SNTP) -RBP* (SNQP*CS221 

if +CSQP*SN221*SNR2D2) -RL*SND2*SN221*CSR2D2+CSPHP* (RCP* 

if CSPHP+RL*SNPHP*SND2) * (SNTP*CS221-CSTP*SNR2D2*SN221) 

ZF2A=RBP* (SNPHP*SNQP-CSPHP*SNTHP*SND2/CSD2) * (CSPHP*SNR2D2* 
if CSGM2/SNTP-CSPHP*CSR2D2*SNGM2*CS221/SNTP+SNR2D2*SNPHP* 

if SNGM2*CS221*CSTP/SNTP+SNPHP*CSR2D2*CSGM2*CSTP/SNTP+SNPHP* 

if SNGM2*SN221) -RBP*SNGH2* (SNQP*SN221-CSQP*CS221*SNR2D2) 

if +RBP*CSR2D2*CSQP*CSGM2 

ZF2=ZF2A+SNGM2*CSPHP* (SNTP*SN221+CSTP*SNR2D2*CS221) * (RCP*CSPHP 
if +RL*SNPHP*SND2) +SNGM2*CSR2D2*CS221* (RCP*CSPHP*SNPHP-RL* 

if SND2*CSPHP**2) +RCP*CSGM2*CSPHP* (CSPHP*CSR2D2*CSTP-SNPHP* 

if SNR2D2) +RL*CSPHP*SND2*CSGM2* (CSPHP*SNR2D2+SNPHP*CSR2D2* 

if CSTP) +RL*SND2* (SNGM2*CSR2D2*CS221-SNR2D2*CSGM2) 

RNXF2=CSPHP*CSGM2*SNTP*SN221+SNPHP*(SNGM2*SNR2D2+CSGM2*CSR2D2* 
if CS221) -CSPHP*CSTP* (SNGM2*CSR2D2-CSGM2*SNR2D2*CS221) 

RNYF2=CSPHP*CS221*SNTP-SN221*(SNPHP*CSR2D2+CSPHP*SNR2D2*CSTP) 
RNZF2=CSPHP*CSTP* (CSGM2*CSR2D2+SNGM2*SNR2D2*CS221) +SNPHP* ( 
if SNGM2*CSR2D2*CS22 1 -CSGM2*SNR2D2) +CSPHP*SNGM2*SN22 1 *SNTP 

RNXH1=-CSPHF*CSGM1*SNTF*SN111+SNPHF*(SNGM1*SNR1D1+CSGM1* 
if CSR1D1*CS111)+CSPHF*CSTF*(SNGM1*CSR1D1-CSGM1*SNR1D1* 

if CSlll) 

RNYH1=CSPHF*SNTF*CS 1 1 1+SNl 1 1 * (SNPHF*CSR 1D1-CSPHF*SNR IDI *CSTF) 
RNZH1=CSPHF*CSTF*(CSGM1*CSR1D1+SNGM1*SNR1D1*CS111)+SNPHF*( 
if CSGM1*SNR1D1-SNGM1*CSR1D1*CS111)+CSPHF*SNGM1*SNTF*SN111 

RNXF1=(CSDLTA*CSGM1**2+SNGM1**2) *RNXHl+SNDLTA*CSGMl*RNYHl+ 
if (1 .DOO-CSDLTA) *CSGM1*SNGM1*RNZH1 

RNYF1=-SNDLTA*CSGM1*RNXH1+CSDLTA*RNYH1+SNDLTA*SNGM1*RNZH1 
RNZFl* (1 .DOO-CSDLTA) *CSGM1*SNGM1*RNXH1-SNDLTA*SNGM1*RNYH1+ 
if (CSDLTA*SNGM1 **2+CSGMl **2) *RNZH1 

F(1)=XF1-XF2 
F(2)=YF1-YF2 
F(3)=ZF1-ZF2 
F(4)=RNXF1-RNXF2 
F(5)=RNYF1-RNYF2 
RETURN 
END 
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_ ^ ^ ***** SUBROUTINE FCN2 ***** 

• • • 

SUBROUTINE FCN2 (X, F,N5 ,PAR5) 

IMPLICIT REALMS (A-H,0-Z) 

INTEGER N5 

REAL*8 X(N5) ,F(N5) ,PAR5(N5) 

COMMON/Al/CNST 

COMMON/A2/DE 1 , DLl , RCP , RL , RCF, PHFDG 

COMMON/ A3/SNPHF , CSPHF , SNR IDI , CSR IDl , SNBP , CSBP , SNR 1D2 , CSR 1D2 
COMMON /A5/SNGM1 , CSGMl , GMAl , SNGM2 , CSGM2 , SND2 , CSD2 , SNDl , CSDl 
COMMON/ A6/SNTP , CSTP , TP , TPDG , SNTP , CSTF , TF , TFDG 
COMMON/ A7 / SNDLTA , CSDLTA , SNPHP , CSPHP , SNR2D2 , CSR2D2 , DLTADG 
COMMON/ A1 2/UF , BF , RBP , UP , RMFl , RMP2 

COMMON/ A18/SNQP , CSQP , QPDG , SNQF , CSQF , QFDG , QP , QF , QPO , QFO 
COMMON/A27/XF1 , YFl , ZFl ,XF2 , YF2 , ZF2 

COMMON/ A34/ FE 1 , SNFE 1 , CSFEl , FE2 , SNFE2 , CSFE2 , FEP , FE221 , FEl 1 1 , 

// FEPDG,FEFDG,F111DG,F221DG 

COMMON/A35/SNTHP , CSTHP , SNTHF, CSTHF , THPDG, THFDG 
QF=QFO+X(3) 

SNQF=DSIN (QF) 

CSQF=DCOS (QF) 

QP=QPO+FEP 

SNQP=DSIN(QP) 

CSQP=DCOS (QP) 

FE111=X(A) 

FE221=X(5) 

SN221=DSIN(FE221) 

CS221=DCOS(FE221) 

SN111=DSIN(FE111) 

CS111=DC0S(FE111) 

sntf=dsin(x(D) 

CSTF»DCOS (X(D) 

SNTP=DSIN(X(2)) 

CSTP=DCOS(X(2)) 

THF=X(1)-QF 
THP=X(2)-QP 
THFDG=THF/CNST 
THPDG=THP/CNST 
SNTHF=DSIN (THF) 

CSTHF=DCOS(THF) 

SNTHP=DSIN (THP) 

CSTHP=DCOS (THP) 

UF1= ( (RCF*CSPHF/SNPHF-RL*SND1-DL1*SNR1D1/CSR1D1) *CSPHF*SNTF+ 

# BF* (-SNPHF*SNQF+CSPHF*SNTHF* (RMFl-SNRlDl) /CSRlDl) -DEI* ( 

# SNPHF-CSPHF*CSTF*SNR1D1/CSR1D1)) /SNTP 
XH1=(CSGM1*CSR1D1*CS111+SNGM1*SNR1D1)*(RCF*CSPHF/SNPHF-UF1* 

# CSPHF-RL*SND1 ) +BF* (CSGMl * (~SNQF*SN 1 1 1-CSQF*CS 1 1 1 * 

# SNR1D1)+SNGM1*CSQF*CSR1D1)+DL1*(SNGM1*CSR1D1-CSGM1* 

# SNR1D1*CS1 1 1) -SNl 1 1*DE1*CSGM1-UF1*SNPHF* (CSGMl* (SNTP* 

# SNl 1 1+CSTF*CS1 1 1*SNR1D1) ~SNGM1*CSTF*CSR1D1) 
YH1=CSR1D1*SN1 1 1* (RCF*CSPHF/SNPHF-UF1*CSPHF-RL*SND1) ~BF* ( 

// -SNQF*CS1 1 1+CSQF*SN1 1 1*SNR1D1) -DLl*SNRlDl*SNl 1 1+DEl* 

# CS 1 1 1+UF1*SNPHF* (SNTF*CS1 1 1~CSTF*SN1 1 1*SNR1D1) 

ZH1A= (CSGM1*SNR1D1-SNGM1*CSR1D1*CS111)*(RCF*CSPHF/SNPHF-UF1* 
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// CSPHF)+RL*SND1*(SNGM1*CSR1D1*CS111-CSGM1*SNR1D1)-BF*(SNGM1 

# * (-SNQF*SN1 1 1-CSQF*CS1 1 l*SNRlDl) -CSGMl*CSQF*CSRlDl) + 

# DL1*(SNGM1*SNR1D1*CS111+CSGM1*CSR1D1)+DE1*SN111*SNGM1 
ZH1=ZH1A+UF1*SNPHF* (SNGMl* (SNTF*SN1 1 1+CSTF*CS1 1 1*SNR1D1) +CSGM1* 

# CSTF*CSR1D1) 

XF1= (CSDLTA*CSGM1**2+SNGM1**2) *XHl+SNDLTA*CSGMl*YHl+ (1 . ODOO- 
// CSDLTA)*CSGM1*SNGM1*ZH1 

YF 1 =-SNDLTA*CSGM 1 *XH 1 +CSDLTA* YH 1+SNDLTA*SNGM 1 *ZH 1 
ZF1= (1 . ODOO-CSDLTA) *CSGM1*SNGM1*XH1-SNDLTA*SNGM1*YH1+ 

# (CSDLTA*SNGM1**2+CSGM1**2) *ZH1 

XF2A=RBP* (CSPHP*SNTHP*SND2/CSD2+SNPHP*SNQP) * (CSPHP*CSR2D2 
// *CSGM2*CS221/SNTP+CSPHP*SNR2D2*SNGM2/SNTP-SNR2D2*SNPHP* 

# CSGM2*CS221*CSTP/SNTP+SNPHP*CSR2D2*SNGM2*CSTP/SNTP-SNPHP* 

// CSGM2*SN221) +CSPHP*CSGM2* (RCP*CSPHP+RL*SNPHP*SND2) * (SNTP* 

# SN221+CSTP*CS221*SNR2D2) 

XF2=XF2A+RBP*CSGM2* (-SNQP*SN221+CSQP*CS221*SNR2D2) +RL*SND2* 

# (CSGM2*CSR2D2*CS221+SNGM2*SNR2D2)-RBP*CSR2D2*CSQP*SNGM2 

# +CSGM2*CSPHP*CSR2D2*CS221* (RCP*SNPHP-RL*CSPHP*SND2) -RCP* 

// CSPHP*SNGM2* (CSPHP*CSR2D2*CSTP-SNPHP*SNR2D2) -RL*CSPHP*SND2 

# *SNGM2* (SNPHP*CSR2D2*CSTP+CSPHP*SNR2D2) 
YF2=CSPHP*CSR2D2*SN221* (RL*CSPHP*SND2-RCP*SNPHP) +RBP* (-SNPHP* 

# SNQP-CSPHP*SNTHP*SND2/CSD2)*(SNPHP*CS221-SNPHP*SNR2D2* 

// SN221*CSTP/SNTP+CSPHP*CSR2D2*SN221/SNTP)-RBP*(-SNQP*CS221 

# +CSQP*SN221*SNR2D2) -RL*SND2*SN221*CSR2D2+CSPHP* (RCP* 

// CSPHP+RL*SNPHP*SND2)* (SNTP*CS221-CSTP*SNR2D2*SN221) 

ZF2A=RBP* (-SNPHP*SNQP-CSPHP*SNTHP*SND2/CSD2) * (CSPHP*SNR2D2* 

# CSGM2/SNTP-CSPHP*CSR2D2*SNGM2*CS221/SNTP+SNR2D2*SNPHP* 

if SNGM2*CS221*CSTP/SNTP+SNPHP*CSR2D2*CSGM2*CSTP/SNTP+SNPHP* 

if SNGM2*SN221)-RBP*SNGM2*(-SNQP*SN221-CSQP*CS221*SNR2D2) 

if +RBP*CSR2D2*CSQP*CSGM2 

ZF2=ZF2A+SNGM2*CSPHP* (SNTP*SN221+CSTP*SNR2D2*CS221) * (RCP*CSPHP 
if +RL*SNPHP*SND2) +SNGM2*CSR2D2*CS221* (RCP*CSPHP*SNPHP-RL* 

if SND2*CSPHP**2) +RCP*CSGM2*CSPHP* (CSPHP*CSR2D2*CSTP-SNPHP* 

if SNR2D2) +RL*CSPHP*SND2*CSGM2* (CSPHP*SNR2D2+SNPHP*CSR2D2* 

if CSTP) +RL*SND2* (SNGM2*CSR2D2*CS221-SNR2D2*CSGM2) 

RNXF2=CSPHP*CSGM2*SNTP*SN221+SNPHP*(SNGM2*SNR2D2+CSGM2*CSR2D2* 
if CS221) -CSPHP*CSTP* (SNGM2*CSR2D2-CSGM2*SNR2D2*CS221) 

RNYF2=CSPHP*CS221*SNTP-SN221*(SNPHP*CSR2D2+CSPHP*SNR2D2*CSTP) 
RNZF2=CSPHP*CSTP* (CSGM2*CSR2D2+SNGM2*SNR2D2*CS221) +SNPHP* ( 
if SNGM2*CSR2D2*CS221-CSGM2*SNR2D2) +CSPHP*SNGM2*SN221*SNTP 

RNXH1=-CSPHF*CSGM1*SNTF*SN111+SNPHF*(SNGM1*SNR1D1+CSGM1* 
if CSR1D1*CS111)+CSPHF*CSTF*(SNGM1*CSR1D1-CSGM1*SNR1D1* 

if CSlll) 

RNYH1«CSPHF*SNTF*CS111+SN111*(SNPHF*CSR1D1-CSPHF*SNR1D1*CSTF) 
RNZH 1 =CSPHF*CSTF* (CSGMl *CSR iDl+SNGMl *SNR IDl *CS 111) +SNPHF* ( 
if CSGM1*SNR1D1-SNGM1*CSR1D1*CS1 1 1) +CSPHF*SNGMl*SNTF*SNl 1 1 

RNXF1= (CSDLTA*CSGM1**2+SNGM1**2) *RNXH1+SNDLTA*CSGM1*RNYH1+ 
if (1 .DOO-CSDLTA) *CSGMl*SNGMl*RNZHl 

RNYF1=-SNDLTA*CSGM1*RNXH1+CSDLTA*RNYH1+SNDLTA*SNGM1*RNZH1 
RNZF1=(1. DOO-CSDLTA) *CSGM1*SNGM1*RNXH1-SNDLTA*SNGM1*RNYH1+ 
if (CSDLTA*SNGM1**2+CSGM1**2) *RNZH1 

F(1)=XF1-XF2 
F(2)=YF1-YF2 
FD)=ZF1-ZF2 
F(4) =RNXF1-RNXF2 
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F(5)=RNYF1-RNYF2 

RETURN 

END 

***** SUBROUTINE FNORM ***** 

SUBROUTINE FNORM (DN IXTH , DN IXFE , DN 1 YTH , DN 1 YFE , DN IZTH , DN 1 2FE) 
IMPLICIT REAL*8(A-H,0-Z) 

COMMON/A3/SNPHF , CSPHF , SNR IDl , CSR IDl , SNBP , CSBP , SNR 1D2 , CSR 1D2 
COMMON/A6/SNTP , CSTP , TP , TPDG , SNTP , CSTF , TF , TFDG 
C0MM0N/A12/UF,BF,RBP,UP,RMF1,RMP2 

COMMON /A34/FE1 , SNFEl , CSFEl , FE2 , SNFE2 , CSFE2 , FEP , FE221 , FEl 1 1 , 

# FEPDG,FEFDG,F111DG,F221DG 
DN1XTH=-CSPHF*SNFE1*CSTF+CSPHF*CSFE1*SNR1D1*SNTF 
DN1XFE=-SNPHF*SNFE1*CSR1D1/RMF1-CSPHF*(CSFE1*SNTF/RMF1+SNFE1* 

# CSTF)-CSPHF*SNR1D1*(-SNFE1*CSTF/RMF1-CSFE1*SNTF) 

DN 1 YTH=CSPHF*CS FE 1 *CSTF+CSPHF*SNFE 1 *SNR ID 1 *SNTF 
DN1YFE=SNPHF*CSFE1*CSR1D1/RMF1+CSPHF*(-SNFE1*SNTF/RMF1+CSFE1* 

# CSTF) -CSPHF*SNR IDl * (CSFE 1 *CSTF/RMF 1 -SNFE 1 *SNTF) 
DN1ZTH=-CSPHF*CSR IDl *SNTF 

DN 1ZFE=-CSPHF*CSR IDl *SNTF 

RETURN 

END 

***** SUBROUTINE PNORM ***** 

SUBROUTINE PNORM (DN2XTH , DN2XFE , DN2YTH , DN2 YFE , DN2ZTH , DN2ZFE) 
IMPLICIT REAL*8(A-H,0-Z) 

COMMON/A6/SNTP , CSTP , TP , TPDG , SNTP , CSTF , TF , TFDG 
COMMON/A7/SNDLTA, CSDLTA, SNPHP , CSPHP , SNR2D2 , CSR2D2 , DLTADG 
COMMON/ A12/UF , BF , RBP , UP , RMFl , RMP2 

COMMON/ A34/FE1 , SNFEl , CSFEl , FE2 , SNFE2 , CSFE2 , FEP , FE221 , FEl 1 1 , 

# FEPDG,FEFDG,F111DG,F221DG 
DN2XTH=CSPHP*(SNFE2*CSTP-CSFE2*SNR2D2*SNTP) 
DN2XFE=-SNFE2*CSR2D2*SNPHP/RMP2+CSPHP*(CSFE2*SNTP/RMP2+SNFE2* 

# CSTP) +CSPHP*SNR2D2* (~SNFE2*CSTP/RMP2-CSFE2*SNTP) 
DN2YTH=CSFE2*CSPHP*CSTP+SNFE2*SNR2D2*CSPHP*SNTP 
DN2YFE=-CSFE2*CSR2D2*SNPHP/RMP2+CSPHP*(-SNFE2*SNTP/RMP2+ 

# CSFE2*CSTP) -CSPHP*SNR2D2* (CSFE2*CSTP/RMP2-SNFE2*SNTP) 
DN2ZTH=-CSPHP*CSR2D2*SNTP 
DN2ZFE=-CSPHP*CSR2D2*SNTP 

RETURN 

END 

***** SUBROUTINE DUP ***** 

SUBROUTINE DUP (UPP , DUPTHP , DUPFEP) 

IMPLICIT REAL*8(A-H,0-Z) 

COMMON/A2/DEl,DLl,RCP,RL,RCF,PHFDG 

COMMON/A5/SNGM1 , CSGMl , GMAl , SNGM2 , CSGM2 , SND2 , CSD2 , SNDl , CSDl 
COMMON/A6/SNTP , CSTP , TP , TPDG , SNTP , CSTF , TF , TFDG 
COMMON/ A7/SNDLTA, CSDLTA, SNPHP , CSPHP , SNR2D2 , CSR2D2 , DLTADG 
COMMON/Al2/UF,BF,RBP,UP,RMFl,RMP2 

COMMON/A18/SNQP , CSQP , QPDG , SNQF , CSQF , QFDG , QP , QF , QPO , QFO 
COMMON/A35/SNTHP , CSTHP , SNTHF , CSTHF , THPDG , THFDG 
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UPP=RBP* (SNPHP*SNQP-CSPHP*SNTHP*SND2/CSD2) /SNTP+CSPHP* (RCP* 

# CSPHP/SNPHP+RL*SND2) 

DUPTHP=RBP*SNQP*(-SNPHP*CSTP+CSPHP*SND2/CSD2) /SNTP**2 
DUPFEP=-RBP*SNTHP* (SNPHP-CSPHP*CSTP*SND2/CSD2) /SNTP**2 
RETURN 
END 

***** SUBROUTINE DUF ***** 

SUBROUT INE DUF (UFF , DUFTHF , DUFFEF) 

IMPLICIT REAL*8(A-H,0-Z) 

COMMON /a 2 /DEI , DLl , RCP , RL , RCF , PHFDG 

COMMON/A3/SNPHF , CSPHF , SNRIDI , CSRIDI , SNBP , CSBP , SNR1D2 , CSR 1D2 
COMMON/A5/SNGM1 , CSGMl , GMAl , SNGM2 , CSGM2 , SND2 , CSD2 , SNDl , CSDl 
COMMON/A6/SNTP , CSTP , TP , TPDG , SNTP , CSTF , TF , TFDG 
COMMON/Al 2/UF , BF , RBP , UP , RMFl , RMP2 

COMMON/Al 8/SNQP , CSQP , QPDG , SNQF , CSQF , QFDG , QP , QF , QPO , QFO 

COMMON/A35/SNTHP , CSTHP , SNTHF , CSTHF , THPDG , THFDG 

UFF= (BF* (SNPHF*SNQF+CSPHF*SNTHF* (RMFl/CSRlDl-SNRlDl/CSRlDl) ) - 

# DEI* (SNPHF-CSPHF*CSTF*SNR1D1/CSR1D1) +CSPHF*SNTF* (RCF* 

# CSPHF/SNPHF-RL*SND1-DL1*SNR1DI/CSR1D1)) /SNTP 

DUFTHF= (-BF*SNQF* (SNPHF-'CSTF+CSPHF* (RMFl/CSRlDl-SNRIDl/CSRlDI 

# ) ) -DEI* (CSPHF*SNR1D1/CSR1DI-SNPHF*CSTF) ) /SNTF**2 
DUFFEF* (SNPHF+CSPHF*CSTF* (RMFl /CSR IDl-SNRlDl /CSR IDI) ) * (-BF* 

# SNTHF/SNTF**2) +DE1* (SNPHF*CSTF-CSPHF*SNR1D1/CSR1D1) / 

# SNTF**2 
RETURN 

END 

***** SUBROUTINE DR2 ***** 

SUBROUTINE DR2 (UPP , DUPTHP , DUPFEP , DX2THP , DX2FEP , DY2THP , DY2FEP , 

# DZ2THP , DZ2FEP) 

IMPLICIT REAL*8(A-H,0-Z) 

COMMON/ A2/DE1 , DLl , RCP , RL , RCF , PHFDG 

COMMON/A5/SNGM1 , CSGMl , GMAl , SNGM2 , CSGM2 , SND2 , CSD2 , SNDl , CSDl 
COMMON/A6/SNTP , CSTP , TP , TPDG , SNTP , CSTF , TF , TFDG 
COMMON/ a 7 /SNDLTA , CSDLTA, SNPHP , CSPHP , SNR2D2 , CSR2D2 , DLTADG 
COMMON/Al 2/UF , BF , RBP , UP , RMFl , RMP2 

COMMON/Al 8/SNQP , CSQP , QPDG , SNQF , CSQF , QFDG , QP , QF , QPO , QFO 
COMMON/A34/FE1 , SNFEl , CSFEl , FE2 , SNFE2 , CSFE2 , FEP , FE221 , FEl 1 1 , 

# FEPDG,FEFDG,F111DG,F221DG 
DX2THP=-CSPHP*CSFE2*CSR2D2*DUPTHP+SNPHP*SNFE2*(UPP*CSTP+SNTP 

# *DUPTHP) +CSFE2*SNR2D2*SNPHP* (-UPP*SNTP+CSTP*DUPTHP) 
DX2FEP=-RCP*CSR2D2*SNFE2* (CSPHP/SNPHP) /RMP2-CSPHP*CSR2D2* ( 

# -UPP*SNFE2/RMP2+CSFE2*DUPFEP)-RL*SND2*CSR2D2*SNFE2/RMP2 

# -RBP* (-SNFE2*CSQP+SNQP*CSFE2/RMP2) +RBP*SNR2D2* (CSFE2* 

# SNQP-CSQP*SNFE2/RMP2) +SNPHP* (UPP* (SNFE2*CSTP+SNTP*CSFE2 

# /RMP2)+SNTP*SNFE2*DUPFEP) 

DX2FEP=DX2FEP+SNPHP*SNR2D2* (UPP* (-CSFE2*SNTP-CSTP*SNFE2/RMP2) + 

// CSFE2*CSTP*DUPFEP) 

DY2THP=SNFE2*CSR2D2*CSPHP*DUPTHP+SNPHP*CSFE2*(SNTP*DUPTHP+UPP 

# *CSTP) -SNPHP*SNFE2*SNR2D2* (-UPP*SNTP+CSTP*DUPTHP) 
DY2FEP=-RCP*CSR2D2*CSFE2* (CSPHP/SNPHP) /RMP2+CSPHP*CSR2D2* (UPP* 

# CSFE2/RMP2+SNFE2*DUPFEP)-RL*SND2*CSR2D2*CSFE2/RMP2-RBP* 
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# (-CSFE2*CSQP-SNFE2*SNQP/RMP2) -RBP*SNR2D2* (SNFE2*SNQP+ 

y/ CSQP*CSFE2/RMP2) +SNPHP* (UPP* (CSFE2*CSTP-SNFE2*SNTP/RMP2 

y/ ) +CSFE2*SNTP*DUPFEP) 

DY2FEP=DY2FEP-SNPHP*SNR2D2* (UPP* (-SNFE2*SNTP+CSTP*CSFE2/RMP2) 
y/ +SNFE2*CSTP*DUPFEP) 

DZ2THP=SNR2D2*CSPHP*DUPTHP+SNPHP*CSR2D2*(-UPP*SNTP+CSTP* 
y/ DUPTHP) 

DZ2FEP=CSPHP''SNR2D2*DUPFEP+RBP*CSR2D2*SNQP+SNPHP*CSR2D2*(-UPP 
y/ *SNTP+CSTP*DUPFEP) 

RETURN 

END 

C.. . 

***** SUBROUTINE DR 1 ***** 

C. . . 

SUBROUTINE DRl (UFF,DUFTHF,DUFFEF,DX1THF,DX1FEF,DY1THF,DY1FEF, 
y/ DZITHF.DZIFEF) 

IMPLICIT REAL*8(A-H,0-Z) 

COMMON/A2/DE1 , DLl , RCP , RL , RCF , PHFDG 

COMMON/A3/SNPHF , CSPHF , SNR IDl , CSR IDl , SNBP , CSBP , SNR1D2 , CSR1D2 
COMMON/A5/SNGM1 , CSGMl , GMAl , SNGM2 , CSGM2 , SND2 , CSD2 , SNDl , CSDl 
COMMON /A6/ SNTP , CSTP , TP , TPDG , SNTP , CSTF , TF , TFDG 
COMMON/A12/UF , BF , RBP , UP , RMFl , RMP2 

COMMON/A18/SNQP , CSQP , QPDG , SNQF , CSQF , QFDG , QP , QF , QPO , QFO 
COMMON/A34/FE1 , SNFEl , CSFEl , FEZ , SNFE2 , CSFE2 , FEP , FE221 , FEl 1 1 , 

# FEPDG,FEFDG,F111DG,F221DG 
DX1THF=-CSPHF*CSFE1*CSR1D1*DUFTHF-SNPHF*SNFE1*(UFF*CSTF+SNTF* 

y/ DUFTHF) -SNPHF*CSFE1*SNR1D1* (-UFF*SNTF+CSTF*DUFTHF) 

DX1FEF=-RCF*CSR1D1*SNFE1* (CSPHF/SNPHF) /RMF1-CSPHF*CSR1D1* ( 
y/ -UFF*SNFE1/RMF1+CSFE1*DUFFEF)+RL*SND1*CSR1D1*SNFE1/RMF1 

# +BF*(-SNFE1*CSQF+SNQF*CSFE1/RMF1)-BF*SNR1D1*(CSFE1*SNQF 

# -CSQF*SNFE 1 /RMF 1 ) +DL 1 *SNFE 1 *SNR IDl /RMF 1 -DE 1 *CSFE 1 /RMF 1 
DX1FEF=DX1FEF-SNPHF* (UFF* (SNFE1*CSTF+SNTF*CSFE1/RMF1) + 

y/ SNFE1*SNTF*DUFFEF) -SNPHF*SNR1D1* (UFF* (-CSFEl*SNTF-CSTF* 

y/ SNFE1/RMF1)+CSFE1*CSTF*DUFFEF) 

DY1THF=-CSPHF*SNFE1*CSR1D1*DUFTHF+SNPHF*CSFE1*(UFF*CSTF+SNTF* 

# DUFTHF) -SNPHF*SNFE1*SNR1D1* (-UFF*SNTF+CSTF*DUFTHF) 
DY1FEF=RCF*CSR1D1*CSFE1* (CSPHF/SNPHF) /RMF1-CSPHF*CSR1D1* ( 

# UFF*CSFE1/RMF1+SNFE1*DUFFEF)-RL*SND1*CSR1D1*CSFE1/RMF1- 

# BF*(-CSFE1*CSQF-SNFE1*SNQF/RMF1)-BF*SNR1D1*(SNFE1*SNQF 

y/ +CSFE1*CSQF/RMF1)-DL1*CSFE1*SNR1D1/RMF’ -DE1*SNFE1/RMF1 

DY1FEF=DY1FEF+SNPHF*(UFF*(CSFE1*CSTF-SNFE1*SNTF/RMF1)+CSFE1* 

# SNTF*DUFFEF) -SNPHF*SNR1D1* (UFF* (-SNFE1*SNTF+CSFE1*CSTF/ 

# RMF1)+SNFE1*CSTF*DUFFEF) 

DZ 1THF=-CSPHF*SNR IDl *DUFTHF+SNPHF*CSR IDl * (-UFF*SNTF+CSTF* 

# DUFTHF) 

DZ1FEF=-CSPHF*SNR1D1*DUFFEF+BF*CSR1D1*SNQF+SNPHF*CSR1D1*( 
y/ -UFF*SNTF+CSTF*DUFFEF) 

RETURN 

END 

C. . . 

^ ^ ***** SUBROUTINE PRINCIPAL CURVATURE ***** 

C. . . 

SUBROUTINE PC (D1 , D2 , D3 , D4 , D5 , D6 , D7 , D8 , XI , X2 , RKl , RK2) 

IMPLICIT REAL*8(A-H,0-Z) 

A“D1*D7-D5*D3 
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B=D1*D8+D2*D7-D5*D4-D6*D3 

C=D2*D8-D6*D4 

Xl= (-B-DSQRT (B**2-4 . 0D00*A*C) ) / (2 . 0D00*A) 

X2= (-B+DSQRT (B**2-4 . 0D00*A*C) ) / (2 . 0D00*A) 

RK1»- (D1*X1+D2) / (D3*X1+D4) 

RK2=- (D1*X2+D2) / (D3*X2+D4) 

RETURN 

END 

* • m 

^ ^ ***** SUBROUTINE PRINCIPAL DIRECTION ***** 

SUBROUTINE PD (Al , A2 , A3 , A4 , A5 , A6 , PCO , PDl , PD2 , PD3) 

IMPLICIT REAL*8(A-H,0-Z) 

PD1=A1*PC0+A2 

PD2=A3*PC0+A4 

PD3=A5*PC0+A6 

RETURN 

END 

_ ^ ^ ***** SUBROUTINE UNIT VECTOR OF PRINCIPAL DIRECTION **'■ 

SUBROUTINE UNIT (El , E2 , E3 , UNITO , UNITl , UNIT2 , UNITS) 

IMPLICIT REAL*8(A-H,0-Z) 

UNITO=DSQRT (El**2+E2**2+E3**2) 

UNIT1=E1/UNIT0 

UNIT2=E2/UNIT0 

UNIT3=E3/UNIT0 

RETURN 

END 

• • • 

^ ^ ^ ***** SUBROUTINE UTRANl ***** 

SUBROUTINE UTRAN 1 (B 1 , B2 , B3 , UT 1 1 , UT2 1 , UT3 1 ) 

IMPLICIT REAL*8(A-H,0-Z) 

C0MM0N/A5/SNGM1 , CSGMl , GMAl , SNGM2 , CSGM2 , SND2 , CSD2 , SNDl , CSDl 

UT11=B1*CSGM1+B3*SNGM1 

UT21»B2 

UT3 1=B1* (-SNGM1) +B3*CSGM1 

RETURN 

END 

• • • 

. _ ***** SUBROUTINE UTRAN2 ***** 

• • • 

SUBROUTINE UTRAN2 (C 1 , C2 , C3 , UT 1 2 , UT22 , UT32) 

IMPLICIT REAL*8(A-H,0-Z) 

COMMON/A5/SNGM1 , CSGMl , GMAl , SNGM2 , CSGM2 , SND2 , CSD2 , SNDl , CSDl 

UT 1 2=C 1 *CSGM2-C3*SNGM2 

UT22=C2 

UT32=C1*SNGM2+C3*CSGM2 

RETURN 

END 

• • • 

^ ^ , ***** SUBROUTINE ELLIPSE AXIS ***** 

• • • 

SUBROUTINE AXIS (SIGMDG, AXISA, AXISB , ALPH, ALPHDG) 
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IMPLICIT REAL*8(A-H,0-Z) 

REAL*8 SIGM(4) 

COMMON/Al/CNST 

COMMON /A36 / RK12 , RK22 , RKl 1 , RK2 1 , UTXl 1 , UTXl 2 , UTY 1 1 , UTY 12 , UTZl 1 , 

# UTZ12 , UTX2 1 , UTX22 , UTY2 1 , UT Y22 , UTZ2 1 , UTZ22 
COMMON/ A37/RNXF2 , RNYF2 , RNZF2 , RNXF 1 , RN YF 1 , RNZ F 1 
COMMON/a38/DEF,SIGMGG,ALPHGG,AXIAA,AXIBB, ERROR 
AK2=RK12+RK22 

AK1=RK11+RK21 
G1=RK11-RK21 
G2=RK12-RK22 
C... FOR LEFT-HAND GEAR 

CS2 1 = (UTX 1 1 *UTX 1 2+UTY 1 1 *UT Y 1 2+UTZ 1 1 *UTZ 1 2) 

IF(II.EQ.2)G0 TO 1208 

SN21=UTX1 1*RNYF2*UTZ12+UTX12*RNZF2*UTY1 1+RNXF2*UTY12*UTZ1 1- 
// (UTX12*RNYF2*UTZ11+RNXF2*UTY11*UTZ12+UTX11*RNZF2*UTY12) 

GO TO 1209 

1208 SN21=-(UTX11*RNYF2*UTZ12+UTX12*RNZF2*UTY11+RNXF2*UTY12*UTZ11) 

// + (UTXl 2*RNYF2*UTZ 1 1+RNXF2*UTY 1 1 *UTZ 1 2+UTXl 1 *RNZF2*UTY 1 2) 

1209 SIGM21=DATAN2(SN21,CS21) 

5 CALL ANGLE (SIGM21) 

CSSIGM=DC0S(SIGM21) 

SNSIGM=DSIN (SIGM21) 

S I GMDG=S I GM2 1 /CNST 

AA= (AK1-AK2-DSQRT(G1**2-2.0DOO*G1*G2''DCOS(2.DOO*SIGM21)+G2**2) 

# )/4.0D00 

BB= (AK1-AK2+DSQRT (Gl**2-2 . DOO*G1*G2*DCOS (2 . D00*SIGM21) +G2**2) ) 

# /4.0D00 

AXISA=DSQRT (DABS (DEF/AA) ) 

AXISB=DSQRT (DABS (DEF/BB) ) 

R AT 1 0= AX I S A/ AX I S B 
FF=G1*DSIN(2.0D00*SIGM21) 

HH=G2-G 1 *DCOS (2 . D00*S IGM2 1 ) 

ALPH2=DATAN2 (FF,HH) 

ALPH=ALPH2/2 . OODOO 
CALL ANGLE (ALPH) 

ALPHDG=ALPH/CNST 

RETURN 

END 


***** SUBROUTINE WRITE2 ***** 

SUBROUTINE WRITE2 
IMPLICIT REAL*8(A-H,0-Z) 

COMMON/A6/SNTP , CSTP , TP , TPDG , SNTF , CSTF , TF , TFDG 
COMMON/A27/XF1 , YFl , ZFl , XF2 , YF2 , ZF2 

COMMON/A34/FE1 , SNFEl , CSFEl , FE2 , SNFE2 , CSFE2 , FEP , FE221 , FEl 1 1 , 

# FEPDG,FEFDG,F111DG,F221DG 

COMMON /A35/SNTHP , CSTHP , SNTHF , CSTHF , THPDG , THFDG 

COMMON / A36/RK1 2 , RK22 , RKl 1 , RK2 1 , UTX 1 1 , UTXl 2 , UT Y 1 1 , UT Y 1 2 , UTZ 1 1 , 

# UTZ 12 , UTX21 , UTX22 , UTY2 1 , UTY22 , UTZ2 1 , UTZ22 
COMMON / A38/DEF , S I GMDG , ALPHDG ,AXISA,AXISB, ERROR 
WRITE (6, 70002) 

70002 FORMAT (IH ,'**RESULT OF KINEMATIC ERROR & BEARING CONTACT**'/) 
WRITE (6 , 70003) FEPDG , FEFDG , THPDG , THFDG , TPDG , TFDG , F22 IDG , F 1 1 IDG 
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70003 FORMAT (IH FEPDG=' ,G20.12, lOX, 'FEFDG=' ,G20. 12, ' 

# IH,’ THPDG=' ,G20.12,10X, 'THFDG=' ,G20.12,/ 

// IH TPDG =' ,G20.12, lOX, 'TFDG =' ,G20, 12,/ 

// IH , ' F221DG=' ,G20.12,9X, 'F111DG=' ,G20.12) 

WRITE (6 , 70004) XFl , YFl , ZFl , XF2 , YF2 , ZF2 

70004 FORMATdH , ' XFl= ' ,G20. 12 , 3X, ' YFl= ' , G20. 12 , 3X, ' ZFl= ' , G20. 12 , / 

// IH , 'XF2=' ,G20.12,3X, ’YF2=' ,G20.12,3X, 'ZF2=' ,G20.12/) 

WRITE (6, 70005) ERROR 

70005 FORMATdH , 'KINEMATIC ERROR' ' ,G20. 12/) 

WRITE (6, 70006) 

70006 FORMATdH , 'RESULT OF PRINCIPAL CURVATURE') 

WRITE (6 , 70007) RKl 2 , RK22 , RKl 1 , RK2 1 

70007 FORMATdH ,'GEAR : RK12= ' ,G20. 12 , 3X, ' RK22= ' , G20. 12 , / 

// IH, 'PINION: RK11=',G20.12,3X, 'RK21=' ,G20.12/) 

WRITE(6,70008)UTX11,UTY11,UTZ11,UTX21,UTY21,UTZ21 

70008 FORMATdH , 'PINION UNIT VECTOR OF PRINCIPAL DIRECTION ',/ 

# IH , 'UTX11=' ,G20.12, 'UTY11=' ,G20. 12, 'UTZ11=' ,G20. 12,/ 

# IH , 'UTX21=' ,G20.12, 'UTY21=' ,G20. 12, 'UTZ21=' ,G20. 12/) 
WRITE (6 , 70009) UTX12 , UTY 12 , UTZ 12 , UTX22 , UTY22 , UTZ22 

70009 FORMATdH ,'GEAR UNIT VECTOR OF PRINCIPAL DIRECTION',/ 

# IH , 'UTX12=' ,G20.12, 'UTY12=' ,G20.12, 'UTZ12=' ,G20.12,/ 

# IH , 'UTX22=' ,G20.12, 'UTY22=' ,G20.12, 'UTZ22=' ,G20.12/) 

WRITE (6, 70010) SIGMDG,ALPHDG,AXISA,AXISB 

70010 FORMATdH , 'DIRECTION & DIMENSION OF ELLIPSE',/ 

# IH , 'SIGMDG=' ,G20. 12, lOX, 'ALPHDG=' ,G20. 12,/ 

# IH , 'AXISA=' ,G20. 12, IIX, 'AXISB=' ,G20. 12,/ 

// IH ,72('.')/) 

RETURN 

END 

C. . . 

***** SUBROUTINE HEADl ***** 

C. . . 

SUBROUTINE HEADl (JJ, II) 

IMPLICIT REAL*8(A-H,0-Z) 

WRITE (6, 60000) 

WRITE (6, 60001) 

WRITE (6, 60002) 

IF(JJ .EQ. DTHEN 
WRITE (6, 60005) 

ELSE 

WRITE (6, 60006) 

END IF 

WRITE (6, 60002) 

IF(II.EQ.2) GO TO 1 
WRITE (6, 60003) 

GO TO 2 

1 WRITE (6, 60004) 

2 WRITE (6,60002) 

WRITE(6, 60001) 

WRITE (6, 60000) 


60000 

FORMAT (IHI,' ') 




60001 

FORMAT (IH * 

: -k k k k k k k k k k k k 

: k k k k k k k k k k k k k k k k k k k k k k k k k k » ^ 

60002 

FORMATdH ['* 



*.) 

60003 

FORMATdH ,'* 

CONVEX 

PART ANALYSIS 

*.) 

60004 

FORMATdH ,'* 

CONCAVE 

PART ANALYSIS 

*.) 
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60005 FORMAT (IH LEFT-HAND GEAR 

60006 FORMAT (IH RIGHT-HAND GEAR 

RETURN 

END 

C. . . 

C • * * 

Q ***** WRITEO(II) ***** 

C. . . 

SUBROUTINE WRITEO(II) 

IMPLICIT REAL*8(A-H,0-Z) 

COMMON/ A2/DE 1 , DL 1 , RCP , RL , RCF , PHFDG 
COMMON/A6/SNTP , CSTP , TP , TPDG , SNTP , CSTF , TF , TFDG 
COMMON/A33/PHP 1 , PHFl , PHP2 , PHF2 , PHP IDG , PHFIDG , PHP2DG , PHF2DG 
IF(II.EQ.2)G0 TO 1 
PP=PHP1DG 
PF=PHF1DG 
GO TO 90000 
1 PP=PHP2DG 
PF=PHF2DG 

90000 WRITE (6, 90001) 

90001 FORMATdH , 'GEAR PARAMETER' ,27X, 'PINION PARAMETER') 

WRITE (6 , 90002) PP , PF , RCP , RCF , TPDG , TFDG 

90002 FORMATdH , ' PHPDG= ' ,G20. 12, 15X, ' PHFDG= ' ,G20. 12, / 

// 1H,'RCP =' ,G20.12,15X, 'RCF “'.G20.12,/ 

// IH ,'TPDG =' ,G20.12,15X, 'TFDG =',G20.12) 

RETURN 

END 
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